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Abstract 

Speciation of mixed lignad complexes of Co(n), Cu(II) and Zn(n) with L-omithine and L- 
glutamic acid has been studied in varying concentrations (0-36.%) of urea-water solutions 
maintaining an ionic strength of 0.16 mol dm~^ (KCl) at 303 K. Titrations were carried out in 
the presence of different relative concentrations of metal (M) to L-ormithine (L) and L- 
glutamic acid (X) (M:L:X = 1:2:2, 1:2:4, 1:4:2) with potassium hydroxide. Stability constants 
of ternary complexes were refined with M1MQUAD75. The best-fit chemical models were 
selected based on statistical parameters and residual analysis. The species detected aieMLXH, 
MLXH^ , MLXH 2 and MLKH^* . The extra stability of ternary complexes compared to their 
binary complexes is believed to be due to electrostatic interactions of the side chains of ligands. 

(Keywords : mixed ligand complexes/L-ormithine/L-glumamic acid/computer modeling/urea). 
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Introduction 

The specificity and selectivity of enzyme-substrate reactions are achieved in vivo 
by manipulating the equivalent solution dielectric contant’ at the active site. The 
equivalent solution dielectric constants of the active sites in bovine carobonic 
anydrase and carboxypeptidase were estimated to be 35 and 70, respectively. In 
carbonic anhydrase the binding of water (a substrate) to the zinc ion simply holds it in 
the correct stereochemical position suitable for attack by CO 2 , where as in 
carboxypeptidase the binding of the substrate results in polarization thus facilitating 
hydrolysis. This variation is brought out by the changes in the interactions of the side 
chains of the protein moiety among themselves and with those of the solvent 
molecules. Further, intramolecular and ligand-ligand stacking interactions in mixed 
ligand complexes are favoured in water-organic media. Thus, the resulting reduced 
equivalent solution dielectric constant, at these active sites, is a consequence of solute- 
solvent and dipole-dipole interactions between the ligands and hence, knowledge of 
the equivalent solution dielectric constant at the active site can throw light on the 
mechanism of the reaction. As a result the protonation constants and binary metal- 
ligand stability constants have been studied^'® in DMF-water and urea-water mixtures. 

The hepatic mitochondrial enzyme, ornithine aminotransferase catalyses the 
transfer of amino group of L-omithine to a-oxoglutarate to form glutamic acid’. The 
enzyme serves an important metabolic fimction in the regulation of ornithine available 
for participation in the urea cycle*’. Hence, in this communication we report the 
stability constants of ternary complexes of Co(II), Cu(II) and Zn(II) with L-ormithine 
and L-glutamic acid, in media of comparable dielectric constant with physiological 
fluids. 


Materials and Method 

Solutions of L-omithine.HCl and L-glutamic acid (E Merck, Germany), Co(II)-, 
Cu(II)- and Zn(II)- chlorides were prepared by dissolving in triple distilled water. 
Urea was recrystallised twice fi-om water and was dried at 60°C for 2 hours'®. Stock 
solutions of urea was prepared by dissolving the appropriate quantity of the purified 
sample in water and stored frozen. The solution was never allowed to stand at room 
temperature continuously for more than 24 hours. 

The titrations were carried out in the medium containing varying amounts of 

urea maintaining an ionic strength of0.16 mol dm"* with KCl at 303. 0±0.1 K. An 
ELICO LI-120 pH meter of 0.01 readability was used. The glass electrode was 
equilibrated in a well-stirred urea-water mixture containing inert electrolyte. 
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Titration of strong acid with alkali was carried out at regular intervals to check 
whether complete equilibrium was achieved. The calomel electrode was refilled 
with urea-water mixture of equivalent composition as that of the titrand. In each of 
the titrations the titrand consisted of mineral acid of approximately 1 mmol in a 
total volxxme of 50 cm^. Titrations were carried out in presence of different relative 
concentrations of metal (M) to L-ornithine (L) and L-glutamic acid (X) (M:L:X = 
1 : 2 ; 2, 1 ; 2 : 4, 1 ; 4 : 2) with KOH. The analytical concentrations of the 
ingredients are given in Table 1 . 

Table 1 -Total initial concentrations of ingredients (in mmol) for mixed-ligand titrations in urea-water 
mixtures. 


[KOH] = 0.4 mol dm"^ Vo = 50.0 cm^ Temp. = 303 K 
Mineral acid =1.00 nmol; p = 0. 1 6 mol dm"^ 


%W/WUrea 

0.00 

Co(n) 

0.143 

TMO 

Cu(n) 

0.147 

Zn(n) 

0.155 

TLO 

L-ORN 

0.307 

0.307 

0.613 

L-GLU 

0.308 

0.616 

0.308 

M:L:X 

1:2:2 

1:2:4 

1:4:2 

5.80 

0.143 

0.147 

0.155 

0.299 

0.300 

1:2:2 





0.299 

0.600 

1:2:4 





0.598 

0.300 

1:4:2 

11.52 

0.149 

0.152 

0.158 

0.301 

0.300 

1:2:2 





0.301 

0.600 

1:2:4 





0.601 

0.300 

1:4:2 

29.64 

0.149 

0.152 

0.158 

0.301 

0.300 

1:2:2 





0.301 

0.600 

1:2:4 





0.600 

0.300 

1:4:2 

36.83 

0.149 

0.152 

0.158 

0.293 

0.301 

1:2:2 





0.293 

0.601 

1:2:4 





0.586 

0.301 

1:4:2 


Vo = Total Volume 

p = Ionic Strength 

TMO = Total metal concentration 

TLO = Total ligand concentration. 





Table 2 -Best fit chemical models of mixed ligand complexes of M(n) with L-omithine and L-glutamic acid in urea-water mixtures. 
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Table 2 Contd.. 
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Result and Discussion 

A preliminary investigation of alkalimetric titrations of mixtures containing 
different mole ratios of L-omithine and L-glutamic acid, in presence of mineral acid 
and inert electrolyte inferred that no condensed species are formed. 

The formation constants^** for acid-base equilibria and for those binary metal 
complexes of L-omithine and L-glutamic acid were fixed in the refinement of mixed 
ligand' stability constants in testing various chemical models using the computer 
program MINIQUAD75”. The best fit model was chosen as that with low standard 
deviation in formation constants and minimum U (corrected for degrees of freedom) 
which was corroborated by other statistical parameters like R etc., given in Table 
2. The stability constants of the complexes are found to follow the order Co(II) < 
Cu{II) > Zn(II). 

Distribution diagrams : 

A perusal of the distribution diagram (Fig. 1) reveals that at very low pH the 
concentration of mixed ligand complexes are less than those of protonated ligands. As 
the pH increased the concentrations of the ternary species increased. In the case of 
copper complexes the concentration of the ternary species are more than even those of 
the ligand species. MUCH and MLXHl showed unusual distribution patterns with 
increase in pH, whereas MLXHl* showed a maximum. The MLXHl* species may be 
assumed to be formed in all these systems according to the equilibrium. 

M(II) + LHl* + XH; = [M(LHJ(XH) p* + 3H * 

In the pH region 3.0-5 .0 ornithine and glutamic acid exist as LHpsad XHz, 
respectively^. These protonated ligands interact with the metal ion to form MLXHl . 
The same species may also be formed due to the interaction between MLX and XHt or 
due to deprotonation of MLXH^* 

M(n) + lh; + xh2 (lh) {XHf + iiP 

MUf-^ +XH2^ MLXHl 
MLXHf^MLXL[+Lr 
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In the same pH region (3. 0-5.0), the existence of AiLW species can be explained 
based on the deprotonation of MLXH* species and also due to interaction of the metal 
ion with ligand species. Even though the deprotonation of MLXH I and interaction of 
metal ligand species lead to the formation of MLXH species, the tendency for the 
formation of h4LH^* species simultaneously in this pH range, results in a plateau in 
this region. Further, increase in MLXH species is noticed with increase in pH. This 
trend is mainly attributed to the absence of formation of MLH species at higher pH 
values where the other two factors dominate. 

M{H) + LH; + XH2 ^M{LH)XotML (XH) + 2Lt 



Fig. 1 “Species distribution diagrams of Co(II)-ORN-GLU ; 

MLXHI*(3)XH2(,4)XH,(,5) LH^ (6)ZJ7(7)X'-(8)M^L 

But at higher pH (5. 0-8.0), MLH^*, which was dominant in binary complexes, 
interacts with XH form of glutamic acid to form MLXH. Deprotonation of MLXHl 
also takes place producing MLffiT. 

MIJf*+XH-=MiLH)X+H" 
much; =MLXH + Lt 



180 


M. SARATCHANDRA BABU et al. 


The CuLXHl and CxxLXH species are formed according to the above equilibrium 
at lower pH region. Hence, the plateau observed for these species in the distribution 
diagram for Co(II) and Zn(II) system is absent in Cu(II) system. 

The concentrations of binary complexes in the presence of ternary species are less 
than 10% in all the systems. Hence, they are not depicted in the distribution diagrams. 

Stability of ternary complexes ; 

The change in stability of ternary complexes as compared to their binary 
analogues was quantified. In one of the approaches’^'’^ the difference in stability (A 
log K) for the two reactions ML with X and M (aq) with X, where X is a secondary 
ligand, is compared with that calculated purely on statistical grounds (Eq. 1). 

A log X = log X- - log log X- - log X" (1) 

Another approach’^' to quantify the stability of ternary complexes was based on 
the disproportionation constant (logX) given by the Eq. 2. 

log X= 2 log x;^- log - log (2) 

The values of A log Xmi and log Xnn for the formation of ternary complexes of 
M (L-omithine) (L-glutamic acid) are recorded in Table 3. Co(II) and Zn(II) form 
octahedral complexes with L-omithine and L-glutamic acid. So on statistical grounds 
the expected A log X value is — 0.4. Cu(II) is assumed to form square planar 
complexes and hence, the predicted statistical value of A log X is — 0.6. The higher 
values observed in the present study infer that the ternary complexes are more stable 
than the binary complexes. 


Table 3- Variation of stability of temaiy complexes of M(n) (ornithine) (glutamic acid) in urea-water 
mixtures. 


%W/WUrea 

Co(n) 


Cu(n) 


Zn(n) 


ALog ATnn 

LogXiin 

A Log Kuii 

LogXim 

ALogiTiiii 

0.00 

8.38 

18.53 

-2.28 

-0.91 

4.82 

5.80 

792 

17.58 

-1.61 

1.36 

6.66 

11.52 

7.31 

16.15 

0.01 

5.87 

5.46 

20.31 

6.80 

15.64 

-1.98 

2.18 

5.37 


Table 3 Contd. 
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Table 3 contd.. 

29.64 6.72 16.01 -0.40 6.66 5.89 

36.83 8.19 17.62 1.34 8.85 6.89 

A log Xino ~ log Piiio - log Pnoo ~ log Pioio 
logXiao = 2 log Piiio - log Pi2oo - log P1020 
A log Klin = log Piiii - log Piioi - log Pioio 
logXiin = 2 log Piiii - log P1202 - log P1020 

The higher values of log X than those expected on the basis of statistical grounds 
(0.6) account for the extra stability of the ternary complexes. Log Xnn for ZnLXH 

species could not be calculated due to the absence of ZnL2-ff ^ species in the Zn(II)- 

Omithine system in urea medium. 

The extra stability of the ternary complexes are due to the interactions outside the 
coordination sphere. This may sometimes be the formation of hydrogen bonds 
between the coordinated ligands. A similar stabilizing effect may likewise be exerted 
by the electrostatic interactions between non-coordinated, charged groups of the 
ligands. Sakurai et carried out extensive studies to establish the laws governing 
interactions of this nature. The higher stability associated with the Metal (L-omithine) 
(L-glutamic acid) species compared to M(L-omithine )2 and M(L-glutamate )2 can be 
explained on similar grounds. It can be attributed to the electrostatic interactions 
between non-coordinated charged groups like-NH*of L-omithine and -COO~ of L- 
glutamate (Fig. 2). The complexation of the metal with two L-amino acids acting as 
bidentate ligands with the resulting electrostatic interaction will result in a MUCH 
species with trans configuration. 



Fig. 1-MLXH species showing electrostatic interactions between the side chains of amino acids. 
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Fig. 3-Variation of log p of Co(II)-ORN-GLU (A), Cu(II)-ORN-GLU (B,C) and Zn(II)-ORN-GLU(D) 
with mole fraction (D-0) and reciprocal of dielectric constant (0-0) in urea-water mixtures. 
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Effect of urea : 

Urea acts as a denaturant of macromolecules (1) by interacting with peptide 
groups through its amido group and (2) by breaking the water structure. The variation 
of stability constants with solvent parameters 1/D or nx is shovwin Fig. 3. As in the 
case of binary complexes, the stabilities of ternary complexes also exhibited non- 
linear trend with increase in urea concentration. For many of the species minima are 
observed at 20.31%. This trend may be due to the complete disruption of water 
structure brought about by urea. The non-linear and irregular variation in the stability 
constants of these ternary complex species may be attributed to complex solute- 
solvent interactions in the aqua-organic mixtures. 

Conclusions 

1. The species confirmed by the present study for complexes of Cu(II) and 
Zn(II) with L-omithine and L-glutamic acid are MUCH, MUCHl and 
MUCH I* . Co(II) formed only MUCH and MUCH]* 

2. The order of the stability of these complexes was determined to be Co(II) < 
Cu(II) > Zn(II). 

3. The concentrations of binary complexes in the presence of ternary species 
are less than 10% in all the system. This and the higher A log K and log X 
values observed in the study infer that the ternary complexes are more stable 
than the binary complexes. The higher stability of ternary complexes can be 
attributed to the statistical factors and electrostatic interactions between non- 
coordinated charged groups like -NH* of L-omithine and -COO“ of L- 
glutamate. 

4. The non-linear trend of stabilities of ternary complexes with increasing urea 
concentration and minima at 20.31% may be due to the complete dismption 
of water stmcture brought about by urea and due to complex solute-solvent 
interactions in the aqua-organic mixtures. 
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Abstract 

Indium(III) is precipitated with oxalic acid in the presence of excess caesium sulphate 
maintaining In^* ; € 204 ^' at 1:2. The chemical analysis of the complex salt obtained 
corresponds to the formula, Cs[In(C 204 ) 2 (H 20 ) 2 ].H 20 . Thermal decomposition studies show 
that the compound decomposes first to the anhydrous caesium indium oxalate and then to the 
final mixture of the oxides through formation of caesium carbonate and indium(III) oxide as 
intermediates. Isothermal study. X-ray diffraction pattern and IR spectral data support the 
proposed thermal decomposition mechanism. 

(Keywords: thermal analysis / bis-oxalates / mdium(in) / caesium salt / X-ray diffraction / IR 
data) 


Introduction 

The thermal decomposition of the oxalates of some transition metal ions and 
lanthanides were investigated by several workers’"®. The presence of suitable 
concentration of oxalates of K^, and Cs"*^ resulted in the precipitation of the 

corresponding complex salts with indium(III) solutions'’"”’. The same precipitates 
were also obtained when indium(III) was titrated directly with oxalic acid in the 
presence of sufficient mono-, di- and tri-valent cations”. Studies on the complex 
formation of the In(ni)-H2C204-H20 system carried out by pH metric”, 
potentiometric’” and solubility methods’^ indicated that the [In(C204)2]' species is 
predominant in the pH range of 2 - 5 . A solid precipitate is, therefore, obvious if 
sufficient precipitating cation is introduced to such a system. The present paper deals 
with the preparation and thermal decomposition studies of the caesium bis- 
oxalatodiaquaindate(III). The complex and its intermediate products of thermal 
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decomposition are characterised using chemical and thermal analyses, infrared 
absorption spectra and X-ray diffraction data. 

Materials and Method 

Thermal analysis unit : SEIKO combined thermal analysis system (TG/DTA-32), 
temperature programmable thermal balance, made in Japan and platinum crucible as 
container were used for taking thermograms in air. The rate of heating is fixed 
tolO°C/inin, and sensitivity of the instrument is 0 . 1 mg. 


Infrared spectra : The infrared spectra of the complexes are recorded on 
SHIMADZU FTIR-8201 PC Infrared Spectrophotometer in KBr pellets. 

X-ray diffraction data : X-ray diffractometer of RICH SEIFERT & CO.(Made in 
Germany) attached to a microprocessor is used for taking X-ray diffraction patterns at 
wave length of Cu Kn = 1 .540598 A. 

The effect of variation of concentration of caesium sulphate was studied on the 
titration of indium(III) with oxalic acid in aqueous medium. It was found that the 
caesium bis-oxalatodiaquaindate(III) complex gets precipitated when the 
concentration of Cs^ is ten times or more than that of In^"^ at the C 20 ffl In^"^ ratio of 
2. Under these conditions the complex was prepared: 

About 50 cm^ of 0.05 M indium(III) sulphate was taken in a 400 cm^ beaker to 
which about 50 cm^ 0.5M caesium sulphate and 100 cm^ triple distilled water were 
added slowly while stirring the contents. Then about 50 cm^, 0.1 M oxalic acid was 
added very slowly (drop wise from a burette) while stirring the contents vigorously. 
The complex salt formed was allowed to settle and filtered through a G 4 sintered glass 
crucible, washed with acidulated water (with H2SO4) to free it from excess of oxalic 
acid, and dried in a vacuum desiccator over silica gel. 


The compound was analysed for its indium, oxalate and water content. 
Indium(ni) was estimated complexometrically using l-(2-pyridyaIazo)-2- 
naphthol(PAN) as indicator''* and oxdate by volumetric titration with standardised 
cerium(IV) sulpahate'^. Water content was determined by difference and from thermal 
data. The results of the analyses are shown in Table 1 . 
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Table l-Cbemicai analysis data of caesium bis-oxalatodiaquaindate(III) 


Methods of 
Preparation 


Composition 

(%) 

Total 

water 

Ratio 

CjO/'/In^ 

Possible Foimula 


In^ 

C204^ 

Cs* 




DeicbmanP] 

— 

— 

— 



Cs[ln(C 204 ) 2 ].nH 20 

Present work 

23.82 

37.63 

27.54* 

11.02* 

2.06 

Cs[In(C204)2].3H20 

Cs[In(C204)2(H20)2].H20 

Calculated 

24.03 

36.85 

27.83 

11.30 

2.00 

Cs[In(C204)2].3H20 


* Calculated from the formula 


Results and Discussion 


Thermal Analysis : 

Thermogravimetric Analysis (IGA) : The thermogram of caesixim bis- 
oxalatodiaquaindate(ni) monohydrate and the data obtained from it are shown in Fig. 1 
and Table 2 respectively. The first loss, which occurs from 30° to 100°C, equals to 
3.77% of the original weight of the complex (observed value is 3.88%) and 
corresponds to the loss of crystal water. This is followed by the complete dehydration 
of the complex salt to give caesium indium oxalate upto ~215°C. The observed weight 
loss is 11.33% against the calcxilated value of 1 1 .30% for the total dehydration of the 
complex. This complex salt eventually exhibits complete dehydration up to around 
200°C. The anhydrous salt remains relatively stable between 200° and 280°C. 


Table 2-Suininary of the thermal decomposition of the caesium salt 


Wt. of 
Comp. 

Step No. 

Temperature 

start °C end °C 

Loss in weight 

Obs. % Calcu. % 

Possible Decomposition product 
(Intermediate) 

12.3 mg 

1. 

30.0 

97.3 

3.88 

3.77 

Cs[In(C204)2(H20)2] 


2. 

150.0 

214.5 

11.33 

11.30 

CsIn(C204)2+ CSpn(C 204 )] 


3. 

214.5 

372.9 

33.87 

33.89 

CS 2 C 2 O 4 + In 203 


4. 

372.9 

418.5 

35.03 

36.82 

CS 2 CO 3 + In 203 





01^5 %/m\n 
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The decomposition of the intennediate(the anhydrous salt) takes place between 
215° and 370°C to give rise to a mixture of caesium oxalate and indium(III) oxide. Of 
course, the actual reaction in this temperature range is none other than the 
decarboxylation of the mdium(III) oxalate which appears to be formed from the 
anhydrous intermediate along with caesium oxalate. The next stage is however, the 
decomposition of the caesium oxalate to the carbonate, which occurs between 370° 
and 420°C. The product at this point is a mixture of CS 2 CO 3 and InaOs that remains 
stable up to nearly 780°C. The above data is in accordance vith the results of the 
earlier study'®. 



5 


-10 


-25 


o 

h- 


*40 


-55 


Fig. 1-TG/DTA curve of caesium bis-oxalatodiaquaindate(I!I) monoliydrate 

Differential Thermogravimetric Analysis (DTG) : Differential losses in weight 
obtained at regular time intervals are noted against temperature for caesium bis- 
oxalatodiaquaindate(in) monohydrate complex as shown in Fig. 1. The losses 
corresponding to the dehydrations of the complex are obvious from the small peaks at 
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approximately 95°C and 170°C. The broad peak with a maximum at 311.0°C, 
covering a range of temperature from 290° to 346°C, corresponds to the 
decomposition of the anhydrous caesium indium oxalate to a mixture of indium(III) 
oxide and caesium oxalate. A small peak around 410°C corresponds to the 
decomposition of caesium oxalate to caesium carbonate. 

Differential Thermal Analysis (DTA) : From Fig. 1 it is clear that the dehydration 
of caesium bis-oxalatodiaquaindate(III) monohydrate takes place in steps. The 
endothermic peak with ATmin at ~90'’C corresponds to the loss of crystal water. A 
relatively strong and broad peak with ATmin at 172.3°C is attributed to the complete 
dehydration of the complex to give the anhydrous caesium indium oxalate. The broad 
exothermic peak is split into three peaks having at 296. 9°C, 313.2°C and 

335. 6°C corresponding to the step-wise decomposition of the caesium indium oxalate 
to give a mixture of caesium oxalate and indium(III) oxide. Another exothermic peak 
with ATmax at 418.4°C signifies the formation of caesium carbonate by decomposition 
of the oxalate. The final product is stable upto 780°C. 

Isothermal decomposition of caesium bis-oxalatodiaquaindate(IlI) monohydrate : 
The thermal decomposition product of the caesium bis-oxalatodiaquaindate(III) 
monohydrate complex salt was obtained by heating ca. 100 mg. of the salt at a 
temperature of 186° C selected from the pyrolysis curve of the complex (Fig.l). The 
decomposition reaction may be described by the equation: 

Cs[In(C204)2(H20)2].H20— ^2:^i£^Cs[ln(C204)2] + 3 H 2 O 

The product so obtained was analysed^’* respectively for indium(III) and oxalate 
contents and the ratio of €204^' to In(III) was found as 2:1. The water content of the 
original complex was also determined as three molecules for each Cs[In(C204)2] 
molecule based on the experimental data. Hence, the possible formula of the caesium 
complex becomes Cs[In(C204)2(H20)2].H20. This is in good agreement with the 
proposed thermal decomposition mechanism of the complex salt given at the end of 
this article. 

Infrared spectra of caesium bis-oxalatodiaquaindate(III) monohydrate 

The infrared spectra of caesium bis-oxalatodiaquaindate(ni) monohydrate and the 
product obtained by heating the complex to 186° C are shown in Fig. 2 and Fig. 3 
respectively. A very strong peak at 3600 cm‘\ and a shoulder at 3477 cm'* and a very 
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broad and strong absorption at 161 1 cm'' (Fig. 2) confirm the presence of water in the 
complex. Similar absorptions, though comparatively less intense in Fig. 3 at 3560 cm'' 
and 1600 cm'' indicate the presence of water in the heated product as well. This water 
in the heated product may be due to absorption of moisture from the atmospheric air. 
A weak absorption at 650 cm'' and a shoulder at 580 cm'' in both figures indicate the 
probable presence of crystal water in the complex as well as in the product obtained 
by heating the complex to 1 86° C. 



Fig. 2-Infrared spectrum of caesium biS“Oxalatodiaquaindate(lII) monohydrate 



A very sharp and strong absorption in the complex at 807 cm'' confirms the 
presence of water in the coorrhnated form. 
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The comparison of the infrared spectral data of caesium bis- 
oxalatodiaquaindate(III) and the heated product with the empirical formula of 
CsIn(C 204)2 are shown in Table 3. 


Table 3- Infrared absorption data of caesium bis-oxalatodiaquaindateflU) monohydrate and the heated 
product obtained at 186°C 



Complex 


Original 

Heated at 186°C 

Band assignment 

3600 sp,vs 

3560 sp,vs 


3470 sh 

3500 sh 

v„.,(H-0-H) 


1680 sh,w 


1611 vb,vs 

1620 b,s 

v.(C=0) + 5(H-0-H) 

1450 m 



1350 sp,m 

1350 sp,m 

Vs(C-O) + v(C-C) 

1300 m 

1300 m 

Vs(C=0) + 5(0-C=0) 

810 sp,s 

810 sp,s 

Coordinated water & 

5 (0-C=0) + V (M-0) 

650 b,w 


Crystal water 

580 sh 



490 b,w 

480 w 

v(M-O) + 5(0-C=0) 

420 sh 

420 sh 


350 vw 

350 w 

5 (0-C=0) + V (C-C) 


Key:- b = broad, m = medium, s = strong, sp = sharp, sh = shoulder, w = weak. 
X-ray Diffraction Data 


The X-ray diffraction data of caesium bis-oxalatodiaquaindate(in) monohydrate 
and that of tihe product obtained after heating the original complex to 186°C and 
cooling are given in Table 4 along with those of indium(III) oxalate for comparison. 

The data in the table clearly show that the d spacings of the heated product does 
not match with the d spacings of In 2 (C 204 ) 3 . So the heated product is not a mixture of 
indium oxalate and caesium oxalate but most probably is a single compound of 
caesium indium oxalate. 
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Table 4 - X-ray diffraction data of caesium bis-oxalatodiaquaindate(III) monohydrate and the product 
after heating it to ISb^C. compared wit h IndiumfUI) oxalate. 

Cs* complex Complex heated to 186 ®C 102(0204)3 


8.3422 



6.8365 


5.6713 

6.3959 

6.41 lx 

5.515x 

6.3I87 


5.357i 

6.1135 


5.0922 

00 

0 


5.0352 

4.8093 


4.794i 

4.5743 

4.6544 

4.648i 

4.5053 


4.0262 

00 

4.2374 

3.92I2 

4.1382 


3.8073 

3.8544 


3.6772 

3.7592 


3.5552 

3.5452 

3.5316 

3.3O61 

3.4982 


3.23I2 

3.47O2 


3.143i 

3.38O5 


3.O5O4 

2.8O62 


3.OOI2 

2.7692 


2.773i 

2.7377 


2.74O2 

2.71 lx 

2.7O85 

2.5292 

2.3692 

2.398i 

2.43I2 

2.3342 


2.3342 

2.3O82 


2.2I81 


Table 4 ConUt.. 
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Table 4 Contd... 


2.2413 


2.I681 

2.1352 


2.134i 

2.0692 


2.0052 

2.0472 

2.0454 

1.952i 

1.9434 

I.95O2 

1.9372 

1.919, 

I.8OI1 

1.840i 

I.74O4 


1.7792 



1.693i 



1.529i 


Further the original complex was also heated to 500°C and maintained at this 
temperature for half an hour. The product obtained in this manner was tested for 
carbonate by the usual acid test. The result of this test indicates the presence of 
carbonate. From the thermal decomposition behavior of the complex it may be 
concluded that caesium carbonate is definitely formed along with the other product, 
indium(III) oxide 

Basing on the results of the above investigations the following thermal 
decomposition mechanism is proposed; 

2Cs[In(C204)2(H20)2]H20 2Cs[In(C204)2(H20)2] 

2 CsIn(C 204)2 ^ CS 2 C 2 O 4 +I 112 O 3 — . > CS 2 CO 3 + In 203 

Overall : 

2Cs[In(C204)2 (H20)2 ].H20 >Cs2C03+In203+6H20+4C0+4C02 

The proposed mechanism suggests the most probable structural formxila of the 
complex with octahedral coordination around indium(III) as : 


Cs[In(C204)2 (H20)2 ].H20. 
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Abstract 

A simple and sensitive visible spectrophotometric method has been described for the assay of 
clindamycin hydrochloride (CM), clomiphene citrate (CP), pimozide (PZ), mefloquine 
hydrochlorde (MQ) or etoposide (EPS) with excess chloromine-T (CAT) and determining the 
consumed CAT with a decrease in color intensity of geillocyanien (GC). All of the variables 
have been optimised. The concentration measurements are reproducible within a relative 
standard deviation of 1 .0%. 

(Keywords : clindamycin hydrocholoiide/clomiphene citrate/pimozide/mefloquine hydrochloride/ 
etoposide/chloramine-T/gallocyanine/spectrophotometiy) 

Introduction 

Clindamycin hydrochloride (CM)'"^, clomiphene citrate (CP)''^ pimozide (PZ)'"^, 
mefloquine hydrochloride and etoposide (EPS)*"^ are useful as antibacterial, 

gonad stimulate, antipsychotic agent, antimalarial and antineoplastic agent 
respectively. The drugs are official in BP* (CM, CP, PZ, MQ and EPS), USP^ (CM, 
CP, PZ and EPS) and IP® (CP). A survey of the literature revealed that only few visible 
(CM’-*°, CP**’'^ PZ*®, MQ*'*’*® and EPS*®’*’) and UV (CP*’® PZ® and MQ®) 
spectrophotometric methods have been reported. The reported spectrophotometric 
method suffer deficiencies such as low s value. It is therefore, of interest to develop a 
simple and sensitive procedure with higher s value for the determination of CM, CP, 



196 


G. P. V. MALLIKARJXJNA RAO et al. 


PZ, MQ, or EPS in pure and pharmaceutical formulations. This paper describes a 
visible spectrophotometric method for the determination of the drug (CM, CP, PZ, 
MQ or EPS) by waldng use of its ability to react with CAT directly. Azine-dyes are 
well known for their high absorptivity and gallocyanine (GC) (phenoxazine-5-ium, 1- 
carboxy-(7-dimethyI-amino)-3, 4-didydroxy chloride : C.I. No. 51030)'* has been 
utilised for estimating excess CAT in the indirect determination of bio-active 
compounds'*’’®. We have applied this sensitive visible spectrophotometric procedure 
in the determination of CM, CP, PZ, MQ or EPS in bulk sample and pharmaceutical 
formulations. The method involves the addition of excess CAT to drug (CM, CP, PZ, 
MQ, or EPS) and unreacted CAT was determined by measurement of the decrease in 
abosorbance of GC. 


Materials and Method 

A Milton Roy Spectronic 1201 UV-Visible spectrophotometer with 1 cm matched 
quartz cells was used for all the absorbance measurements. 

All chemicals were of analytical grade and all solutions were prepared in triply 
distilled water. Aqueous solution of CAT (Loba, 200 pg/ml, 7.1 x lO^M), GC 
(Croma, 100 |xg/ml, 2.969 x lO^M) and hydrochloric acid (E. Merck, 5.0 M) were 
prepared by dissolving the required amounts in triple distilled water. 

Stock solutions of CM, CP, PZ, MQ or EPS (400 |xg/ml) were prepared by 
dissolving 40 mg of drag (CM, CP, PZ, MQ or EPS) initially in either 10 ml of 
distilled water or 0.2 M NaOH (EPS), followed by dilution to 100 ml with distilled 
water. Stock solutions were further diluted stepwise with distilled water to give the 
working standard solutions. 

To each of 25 ml calibrated volumetric flasks containing standard drug solution 
(CM : 0.5-2.5 ml, 80 pg/ml; CP : 0.5-3.0 ml, 100 pg/ml; PZ : 0.5-3.0 ml, 40 pg/ml; 
MQ : 0.5-3.0 ml, 100 pg/ml; EPS : 0.5-3.0 ml, 100 pg/ml), 1.25 ml of 5 M HCl and 2 
ml CAT (200 pg/ml) were added and the solutions were diluted to 20ml with distlled 
water. After 10 min, 5 ml of GC solution (100 pg/ml) was added and mixed 
thoroughly and the absorbances were measured after 5 min at 540 nm against distilled 
water. The blank (omitting drag) and dye (omitting drag and CAT) solution were 
prepared in similar manner and their absorbances were measured against distilled 
water. The decrease in absorbance corresponding to consumed CAT and in turn to 
drag content was obtained by subtracting the decrease in absorbance of the test 
solution (dye-test) from that of the blank solution (dye-blank). The amount of drug in 
a sample was obtained from tire Beer Lambert's plot of drug concerned. 
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An accurately weighed quantity of the drug formulation in the form of capsule 
(CM or EPS), tablet (CP, PZ or MQ) or a measure voliune of injection fluid (EPS) 
equivalent to 40 mg were separately taken fro each drug and extracted with warm 
chloroform (3 x 25 ml) and filtered. The volume of the combined chloroform extract 
was brought to 100 ml with chloroform. The chloform extract of each drug was 
evaporated to dryness and the residue obtained was dissolved to prepare working 
sample solutions, as described in earlier and analysed as described under the 
procedure for pure samples. 


Result and Discussion 

The optimum conditions for the development of the method for each drug were 
established by varying parameters one at a time^° and observing the effect produced 
on the absorbance of the colored species. 

The method involves two steps, namely reaction of drug, (CM, CP, PZ, MQ or 
EPS) with an excess of CAT and the estimation of unreacted CAT using a known 
excess of GC. The excess GC remaining was then measured with a spectrophotometer 
at 540 nm. The effect of CAT concentration and acidity for different time intervals in 
first step and dye concentration in second step, waiting period in each step with 
respect to maximum sensitivity, adhrence to Beer's law, reproducibility and stability 
of final color were studied through control experiments. The observed optimum 
conditions are incorporated in the procedure recommended, the mole ratio studies 
revealed that each mole of the drug consumed 1.8 (CM or CP), 1.37 (MQ), 2.67 (PZ) 
or 1.25 (EPS) moles of CAT. The consumption of CAT by each drug under 
established experimental condition giving reproducible values depends upon the 
nature and impact of functional groups present in it. 

Analytical data : 

The optical characteristics such as Beer's law limits, molar absorptivity and 
Sandell's sensitive for these methods are given in Table 1. The precision of the 
method was found by measuring absorbances of six replicate samples containing 
known amounts of drug and the results obtained are incorporated in Table 1. 
Regression analysis using the method of least square was made to evaluate the 
slope(b), intercept (a) and correlation coefficient® for each drug (Table 1). The 
relative standard deviation and % range of error at 95% confidence level are also 
given in Table 1. The accuracy of the method for each drug was ascertained by 
comparing the result by proposed method and reference method (UV) statistically^° by 
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t- and F-tests (Table 2). As an additional check of accuracy of the proposed method, 
recovery experiments were performed by adding a fixed amount of the drug to the pre 
analysed formulation and the results are presented in Table 2. This comparison show 
that there is no significant difference between the results obtained by the proposed 
method and those of reference. The similarly of the results is an obvious evidence that 
during the application of this method, the excipients that are usually present in 
pharmaceutical formulations do not interfere with the assay by the proposed method. 
Hence, the proposed method is simple, sensitive and useful for the assay of CM, CP, 
PZ, MQ or EPS in pure pharmaceutical formulations. 

Table 1 -Optical characteristics and precision 


Para meter 

CM 

CP 

PZ 

MQ 

EPS 

Wnm) 

540 

540 

540 

540 

540' 

Beefs law limits (pgml“^) 

0.8-8.0 

M2 

0.64.8 

1-12 

2-12 

SandelFs sensitive (pg.cm” 
Vo.OOl absorbance unit) 

0.0157 

0.202 

0.0105 

0.118 

0.0305 

Molar extinction coefficient 
(1. Moie"^ cm-^) 

2.930 X 10“ 

2.960 X 10 

i“ 4.384 X 10“ 2.239X lO’ 

1.92 X 10“ 

Relative standard deviation 

(%)*♦ 

0.738 

0.509 

0.828 

0.585 

0.715 

% range of error 
(confidence limits 0.05 
level) 

0.775 

0.535 

0.869 

0.614 

0.748 

Correlation coefficient <D 

Regression equation Y*: 

0.9999 

0.9999 

0.9999 

0.9999 

0.9999 

Slope (b) 6.325 X 

10"’ 4.94 X 10"’ 

9.5 X 10"’ 

5.331 X 10"’ 

3.30 X 10“^ 

Intercept (a) 2.4 x 10"’ 7.33 x 10"’ 

-8.66 X 10"“ 

1.466 X 10"’ 

-1.067 X 10“^ 


* y = a 4- be where c is concentration 
** six replicate samples 
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Table 2 - Analysis of pharmaceutical formulations by proposed and reference methods 


Drug/Pharmaceutical 

formulations* 

Labelled amount 
(mg) 

Amount foimd** 

Proposed Reference method 

method CM^ CP*•^ PZ^ 

MQ''. EPS* 

% Recovery by 
proposed 
method*** 

CM 


149.75±0.42 



Capsules-I 

150 

F=1.22 

149.710.46 

98.8310.28 



t=l.l 





149.91±0.61 



Capsules-n 

150 

F = 3.17 

148.8911.09 

99.9410.41 



t = 0.55 





148.8410.58 



Capsules-IH 

150 

li 

to 

149.4210.57 

99.2010.38 

CP 


49.9410.50 



Tablets-I 

50 

F= 2.44 

49.8110.32 

99.8811.01 



t = 0.45 





50.0210.25 



Tablets-I 

50 

F = 4.34 

50.1310.12 

100.4410.50 



t = 0.55 





24.9710.097 



Tablets-in 

25 

F = 3.44 





t=l 

25.0510.18 

99.9010.38 

PZ 


1.9810.114 



Tablets-I 

2 

F = 3.72 

2.0010.21 

99.4710.86 



t=1.5 





2.00±0.015 



Tablets-I 

2 

F = 3.07 

2.0111010 

100.2110.37 



t = 2.05 




Table 2 Contd.. 
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Table 2 Contd.. 





2.0010.015 



MQ 

Tablets-in 

2 

F = 2.77 

t= 1.5 

250.6410.59 

1.99±0.025 

99.9110.77 


Tablets-! 

250 

F = 2.27 

t=1.86 

248.9611.17 

250.9710.89 

100.210.23 


Tablets-n 

250 

F = 3.14 

t=1.96 

249.41i:0.59 

250.3810.66 

99.5810.47 

EPS 

Tablets-in 

250 

F= 1.171 

t=1.31 

496.65±1.6 

250.3811.01 

98.8210.225 


Capsules-I 

500 

F=1.03 

t=1.77 

496.6±1.35 

495.111.50 

99.331032 


Capsules-n 

500 

F = 2.27 

t=1.44 

20.19±0.052 

499.112.05 

99.3210.27 


Injection-I 

20 

F = 3.77 

t = 0.21 

20.2010.1 

100.9810.26 


♦ Three types of capsules/tablets (CM, CP, P, MQ or EPS) from three different pharmaceutical 
companies. 


** Average + standard deviation of six determinations; the t~and F- test values refer to comparison 
of the proposed methods with the reference method. Theoretical values at 95% confidence 
limits, t = 2.57, F = 5.05. 

*** Recovery of 10 mg added to the phannaceutical formulations (average of three determinations) 
A Developed in laboratory for CM in O.IN H2SO4 (X^nax 218 nm). 

# Developed in laboratory for EPS in methanol (X^ax 283 nm). 
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Abstract 

Ternary complex formed by Cu(II)-sulphosalicylic acid-glycine system has been studied by 
Cu(n) ion-selective electrode. The value of stability constant (log Kmla = 8.63) determined 
from the measurement of free ion concentration using the Cu(n) ion selective electrode agrees 
with, determined by pH-titration technique (log Kmla = 8.24). 

Introduction 

Ion-selective electrodes have been extensively used for the determination of 
stability constants’. In this laboratory a number of electrodes have been prepared for 
the determination of stability constants of metal complex (1:1) in solution^'^. In this 
paper we describe our attempt to use, for the first time, the ion-selective electrode for 
determination of stability constants of ternary complex (1:1:1) formed by Cu(II)- 
sulphosalicylic acid (SSA)-glycine system. For this a precipitate based Cu(II) ion- 
selective electrode was prepeured as per method given in the literature’. In order to 
compare the values obtained by ion-selective electrode, the system was also studied 
by pH-titration technique. 


Materials and Method 

For that preparation of electrode CuS was used as an electroactive material which 
was obtained by precipitating an acidic CuCl 2 solution by passing H 2 S gas. After 
separating the precipitate from the solution it was washed thoroughly with water and 
dried at room temperature. 100 mg of dried precipitate with an equal amount of 
Araldite (Ciba, India) was homogenously mixed on a piece of filter paper. The paste 
was spread uniformly and allowed to dry in air for 24 hr to form a master disc of ~ 0.5 
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mm thickness. The adhering filter paper was peeled off. A disc of 1 cm diameter was 
cut from the master disc. This was kept immersed in 1.0 Mol dm'^ CuCl 2 solution for 
two days and any portion of the filter paper still adhering to the surface of the disc, 
was removed. The disc dried in air, served as a master membrane. The disc was fixed 
to one end of a glass tube (length, 5 cm) with Araldite and dried. The tube was filled 
with 0.1 mol dm"^ CuCU solution and kept immersed in a solution of 0.01 mol dm"^ 
CuCl 2 for a week. An Ag-AgCl electrode was inserted through the open end of the 
tube for electrical contact. 


A P hili ps H-meter (Model PR 9405 M) was used for potential measurements with 
a calomel electrode as an external reference electrode. Before measurements, the 
prepared membrane electrode was pretreated with a solution of Cu(II) ions. All the 
measurements were carried out at 30 ± T. Double distilled water was used and 
chemicals were of reagent grade. 

The electrode assembly for the measurement of Cu-ion concentration could be 
expressed as 


Ag,AgCl 

CUCI 2 O.OI 

Membrane 

Sample 

External reference 


mol dm"^ 


solution 

electrode 


For the determination of stability constant of binary as well as of ternary 
complexes a titration assembly was made as shown in Fig. 1. The following sets of 
solutions were prepared. 


External reference 
electrode 
(ColomeO 

0-1 mol dni'^ 
Cu CI 2 Soln. 


Membrane 
Mognetrc stirrer 



( Ag-AgQ) Interndl reference 
electrode 


Combined 
Glass electrode 
( pH-meter) 


Beaker 


Test solution 
for titration 


Fig. 1-A schematic representation of titration assembly. 
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1. Metal ion [Cu(II); 0.1 mol dm 2 ml] + primary ligand (sulphosalicylic 
acid; 0.1 mol dm“^, 2 ml) + water. 

2. Metal ion [Cu(II); 0.1 mol dm~^, 2 ml] + secondary ligand (glycine; 0.1 
mol dm"^, 2 ml) + water. 

3. Metal ion [Cu(II); 0.1 mol dm"^, 2 ml] + primary ligand (sulphosalicylic 
acid; 0.1 mol dm”^, 2 ml) + secondary ligand (glycine; 0.1 mol dm”^, 2 ml) 
+ water. 

Since Cu(II) solution also contains an equivalent concentration of the free 
perchloric acid, hence in the system 1,2 and 3, perchloric acid was not added 
separately. Thus the concetration of free perchloric acid was kept constant in all cases. 

In each case the total volume was 50 ml. 

In all these sets the metal ion and ligand concentration was kept constant and pH 
was varied in the range of 3-6.5 by titrating it with 0.2 mol dm’^ of NaOH soltion. The 
pH and the potential of the cell were simultaneously recorded using Century CP 901 
pH meter and Philips pH-meter (Model PR 9405 M) respectively at room temperature 
(30 ± 1°). In both the measurements a calomel electrode was used as a reference 
electrode. 


Results and Discussion 

The study of the mixed ligand complex lies in the potentiometric titration of the 
reaction mixture containing a T.l:l molar ratio of the metal, the primary ligand and 
secondary ligand. The metal and primary ligand as well as metal and secondary ligand 
containing 1:1 molar ratio of metal to ligand are also titrated pH-metrically against 
same alkali solution under identical conditions. 

The titration curve (Fig. 2) for the binary 1:1 Cu(II)-SSA system, a well marked 
inflextion at a « 2. 0-2. 5 indicating the formation of the stable [Cu-SSA] chelate. It 
of interest to note that another well marked inflextion appeared at a « 3. It showed 
the addition of extra OH and possibly formation of a dimeric species M 2 A 2 OH. It may 
also be noted that beyond a « 3 ^H ~7) turbidity appears. In 1:1 Cu(n)-glycine 
system (MA), the titration curve gave the first inflextion at a « 1, followed by a 
second inflextion at a » 2. In this system also the turbidity appeared beyond a » 2 (pH 
« 6.5). It thus shows the formation of ML and ML(OH) respectively. 
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Fig. 2-Titration curves for Cu(n)-Sulphosalicylic-acid-glycine system. 

H: 0.1 MHCIO4. A ;0.1 M glycine, L ; 0.1 M Sulphosalicylic acid, Cu(n) : 0.1 M.; ML: 1:1 molar ratio 
of Cu(II) to sulphosalicylic acid.; MA: 1:1 molar ratio of Cu(n) to glycine.; MLA : 1:1:1 molar 
ratio of Cu(II) to sulphosalicylic acid and glycine. ; C : Composite curve 
; indicates precipitation. 

In 1;1;1 Cu-SSA-glycine mixed ligand system (MLA), the first inflextion 
occurred at a « 4 and the solution remained clear throughout the titration. The 
important observations are also presented in Table 1. 

Moreover the experimental mixed ligand titration curve is well displaced from the 
theoretical composite curve. These observations indicate the formation of the mixed 
ligand complex in the region a = 3-4. 



determination of stability constant of ternary-complex 


207 


Table I- Important observation during pH titrations 



System (titration curve) 

Inflextion points 

Color change and pH 

pH of appearance of 
turbidity 

ML 

a~0-2 

a --2-3 

Very light blue- 

leafy green (pH - 2.9) 
color intensity 
increases. 

pH = 6.8 

MA 

a -0-1 

Light blue color 
intensify 

no turdidity 

MLA 

a -2-3 

Very light blue- 

sea green (pH - 2.6)- 

leaf green (pH - 3.6)- 

bluish green (pH - 6.0)- 

biue 

pH = 6.8 


a = equivalents of alkali added 

The reaction in the simple 1 : 1 Cu-SSA(ML) system may be represented as ; 
Cu^^ + HjL- -» [CuLr+2irat(a = 2-2.5) 

[CuLr + OH--> [CuL(OH)f-at(a = 3) 

The hydroxochelate may further disproportionates as 
[CuL(OH)]^--> [CuL]- + Cu(OH )2 
In the simple 1:1:1 Cu-SSA-glycine system (MLA) 

+ H 2 L — [ML]" + 2¥t at (a = 2-3) 

between 


[ML]- +HA — [MLA]^- + iT (a = 3-4) 
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where represents ion, 

HjL represents dissociated sulphosalicylic acid. 

SSA represents undissociated sulphosalicylic acid. 

HA represents undissociated glycine molecule. 

In the region (a = 3-4), glycine combines with [Cu-SSA]" chelate to form the 
mixed ligand complex [Cu-SSA-gly]^". The result show that in 1:1:1 system, finally 
the mixed complex [MLA]'" is formed. The calculations for stability constants also 
support this contention. The system giving stability constant (Kmla) values which 
agree will within experimental limits. 

The protonation constant of glycine was determined by Iruing and Rossotti's 
method'” in aqueous medium log K” = 9.14 at SO^C, which are in good agreement 
with those reported by Charles and Preiser” 9.78, whereas the dissociation constants 
of sulphosalicylic acid are taken from the previous work of one of the authors'^ as log 

=2.93 and log K” = 12.60 at 30“C. 

In the present study for the potentiometric determination of stability constant of 
1:1:1 mixed ligand system was calculated in terms of the secondary ligand A by the 
methods developed by Martell et a/.'”’'"* and Ozer”, adopted by Singh and Srivastava'® 
which is a modified form of Irving and Rossotti's method. 

The binary complex formation can be shown as 


M+HJL—ML + nH 

(charges are omitted for simplicity) 

The conditional stability constant K\ could be expressed as 

r = M1M. 

‘ mWA 


( 1 ) 


( 2 ) 


However, the stability constant is defined as 
K 


( 3 ) 
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where, 


[ML] K[ 

m\HA (AT 


and (p 


I m I m 
Kd, KdJCd^ 


Kd^Kd^Kd, 



(4) 


Kd^Kd^Kd^ are first, second and third dissociation constants of ligand. 

In case of sulphosaliclylic acid there are two acid groups which dissociate on 
complex formation and these are hydroxyl and carboxyl groups. Therefore, cp can be 
expressed as ; 


(p = 


\ Ag Ag A3 


(5) 


For 1 : 1 complex we can write 

where M° and L° are the total metal and ligand concentrations respectively. 

The equation (6) can be rearranged as 

W] 1 , \M\ (7) 

K[L L 

Substituting K[ from equation (3) in the equation (7) we obtain 

[M] 1_ 1 1 

. M--[M](p K,L cpL- 

If the total metal and ligand concentrations are kept constant and pH of the 
solution is varied, a plot of [M\IK^-[M\. l/cp against [Mj/cp would give a straight line 
with slope equal to 1/L° and intercept equal to \IK\I^. Thus one can calculate the value 
of stability constant of 1 ; 1 complex, Ki from the slope and intercept of such plot. 
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The idea could be extended fiuther and expression for calculation of stability 
constant of 1:1:1 complex could be obtained However, when 1:1:1 complex 
formation starts, the formation of 1:1 complex should be complete. Considered A as a 
secondary ligand, the formation of 1 : 1 : 1 complex can be shown to take place as 

ML + HJ^MLA +nH (9) 

The equilibrium constant for this reaction can be expressed as 


[MLA][H] 

\ML\\HA\ 


The stability constant for this ternary complex can be expressed as 


Where q) has same meaning as defined in equation (4) and (5) but this time the 
dissociation of secondary ligand A is taken into consideration. 

Now the total metal ion concentration can be expressed as 

iW" = [ML] + [MLA] + [Af] (12) 


Therefore, for 1 : 1 : 1 complex, we can modify the equation (8) as 

[M] 11 [M] 1 

{{M- - [ML]-[M]y K^^L-^ (|) L 

Thus a plot of 


[M] 1 . ,\M] 

will give slope as 1/L and intercept as \IK\^\L , from which ATi i can be calculated 
(Fig. 3). 
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Fig. 3~Deteiinination of stability constant of Cu(n)-SSA-glycine complex. 

Plot of M 1 concn. of Cu(n) = 4x10"^ mol. dm~^,;fmal concn. of SSA= 4 x 
PhmFM'cp <p 

10”^ mol. dm'^; final concn. of glycine = 4 x 10”^ mol. dm■^; pH varying, Kdi for glycine = 
7.24 X 10"*°,; temp. = 30 “C, slope = 5.08 xlO"^,; Intercept = 1.18 x 10"’, log Ki = 8.63. 

The electrode gave a linear response down to a cone. 1 x lO"^ mol dm"^. The slope 
of the electrode was 30 mV per decade change in concentration and thi« remains 
constant for over a period of three weeks. 

To find the response time, the electrode was dipped in 0.01 mol dm"^ CuCU 
solution and suddenly it was taken out and dipped in 0.001 mol dm"^ CuCl 2 solution. 
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The values of the potential change were noted every 5s. A constant potential was 
obtained after 30 s. 

To study the effect of pH, the pH of a 0.01 mol dm"^ CuCb solution was varied by 
the addition of NaOH or HCl. It was found that the potential remains unchanged 
within the pH range 3-7.5. Thus the working pH range of the electrode is 3-7.5. 

The pH titration curves given in Fig. 2 reveal that sulphosalicylic acid starts 
forming complex much earlier than glycine and thus in the mixed complex, 
sulphosalicylic acid is a primary ligand and glycine is the secondary ligand. 

Determination of stability constant of ternary complex using ion-selective 
electrode method : After obtaining the free metal concentration from the caliberation 
curve in the pH range 3. 5-6.5, for each pH-value, the graph was plotted using the 
equation (12). The slope and intercept of this plot were determined by the method of 
least squares, which resulted in the value of stability constant (log Kmla) values of 
1:1:1 complex as 8.63. 

Determination of stability constant by pH-titration technique : For pH titration, 
stability constant of ternary complex was also determined using Martell'^’’'* and Ozer'^ 
method (log Kmla = 8.24). 

Thus we see that the value of stability constant of 1:1:1 complex of the ternary 
system Cu(II)- sulphosalicylic acid-glycine determined by ion-selective electrode 
(ISE) agrees quit well with the value of stability constant determined by pH titration 
technique. Therefore, the ion-selective electrode can be used for studying the systems 
where more than one ligand is attached to the metal ion. 

References 

1. Moddy, G.J. & Thomas, J.D.R. (1978) Ion-Selective Electrodes in Analytical Chemistry, Ed. 

Preiser, H. Plenum Press New York, Vol. 1, p. 340 and references therein. 

2. Lai, U.S., Chattopadhyaya, M.C. & Dey, A.K. (1980) J. Indian Chem. Soc., 19(A) : 390. 

3. Lai, U.S., Chattopadhyaya, M.C & Dey A.K. (19SI) Electrochim. Acta 26 : 283. 

4. Lai, U. S. (1982) Doctoral Thesis, University of Allahabad. 

5. Lai, U. S., Chattopadhyaya, M.C & Dey A.K. (1982) J. Indian Chem. Soc. 59 : 493. 

6. Bhattacharya, M. & Chattopadhyaya, M.C ( 1 987) J. Indian Chem. Soc. 64 : 575 . 

7. Bhattacharya, M. (1992) Doctoral Thesis, University of Allahabad. 

8. Kar, R. & Chattopadhyaya, M.C (1991) J. Indian Chem. Soc. 68 ; 459. 

9. Kar, R, Azam, N. & Chattopadhyaya, M.C (1 992) Bull Chem. Soc., Ethiopia 6 : 1 09. 



determination of stability constant of ternary-complex 


213 


10. Irving, H. & Rossotti, H. S. (1953) J. Chem. Soc. 3397 : (1954) 2433. 

11. Charles, R. G. & Frieser, H. (1952) J. Am. Chem. Soc. 14: 1385. 

12. ChattQpadhyaya, M.C (1974) Doctoral Thesis, HT, Bombay. 

13. Carey, G. H., Boguchki, R. F. & Martell, A. E. (1964) Inorg. Chem. 3 : 1288. 

14. L. Heureux, G. A. & Martell, A. E. (1966) J. Inorg. Nucl. Chem. 28 : 481. 

15. Ozer, O. Y. (1970) J. Inorg. Nucl. Chem. 32 : 1279. 

16. Singh, M. K. & Srivastava, M. N. (1973)7. Inorg. Nucl. Chem. 35 ; 2433. 




PROC. NAT. ACAD. SCI. INDIA, Vol. 73, A(n), 2003 


The Magic Square 

JAMUNA PRASAD AMBASHT and STACEY FRANKLIN JONES* 

Mathematics and Computer Science Department Benedict College, 34 Westpine 
Court, Columbia, SC 29212 (U.S.A) 

*Computer Science Department Johns Hopkins University Baltimore, MD 21228 
Mathematics artd Computer Science Department Benedict College, Columbia, SC 29204 
Received July 10, 2001; Accepted January 30, 2002 


Abstract 

The primary purpose is to study the magic square algorithmically, but alongside passing 
reference to magic triangle or such relevant ideas is not ruled out. Chinese claim that the magic 
square was known to them as far back as 2200 B.C. The magic square was known in 
prehistoric days in India as found in stone carvings and inscriptions on walls and gates of old 
ruins. Also found in prehistoric era are some good-luck charms in business records and in 
horoscopes in India. Examination and construction of the magic square provides an 
opportunity for an aspiring mathematician and computer scientist to explore pattern generation 
and respective algorithm development. The magic square is also fertile ground for exploring 
simplification of complex arithmetic manipulation by way of algorithmic characterization. 

1. Introduction 


"I-King" is a Chinese classic. A numerical diagram, 'the Lo-Shu', appears in it. 
Emperor Yu saw, as he walked along the bank of the yellow river, a divine turtle 
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This is a 3X3 magic square 



whose sum is 15, if we add rows or columns. Interestingly enough the sums along the 
leading diagonals are also found to be 15. Thus a systematic study of magic squares 
of odd orders was undertaken. Some rules were developed. 

1 . The rules of arrows: 

a. Always start from the center of the top row. 

b. Extend an arrow diagonally from the top right comer of the cell and 
position the successive integer if possible^^ 

c. If not, allow the arrow-head to reposition to fall to the right and straight 
down'^ inside the square and position the successive integer if possible 

d. Drag the arrowhead all the way to the left come^\ in case (b.) is not 
possible; and the resulting cell at the beginning of the row was still inside 
the square. 

e. If the arrowhead could not follow (b.) or (c.), or whenever the arrow’s 
movement was blocked, position the successive integer in the cell directly 
below the starting position. 

2. The rule of bulging out: 

a. Position dotted cells in a checkered pattern on each of the four sides of the 
square. This creates bulges. 

b. Beginning at the left most point for each diagonal, start imaging 
successive integers (outside the square moving along each diagonal.) 

c. Erase the outside bulges by throwing the number inside far on the opposite 
side. 
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n. We construct a 5X5 magic square. We fix the position of 1 by using 1(a) and 
following the rule of arrows. This is a 5X5 normal magic square. 
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We have bulging out along the dotted lines. We follow rule outlined in 2. 

III. These rules (1. and 2.) were for odd magic squares. It was found that even 
magic squares were much harder and sometimes impossible to construct. 

Now we consider a 4X4 magic square. 


16 

2 

3 

13 

5 

11 

10 

8 

9 

7 

6 

12 

4 

14 

15 

1 


Step 1. Cross out the leading diagonals. 


FINAL ENTRY 



ENTRY A 
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Step 2. Then fill the rest. 

Start filling from the bottom row to top left to right, in reverse order, and continue 
in this fashion with entries that were omitted. 
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Step 1. Cross diagonals of 4 parts of the magic square. Each part is a 4X4 grid. 
Start filling rows leaving blank the places where the diagonals travel. Record your 
blanks in a row. Superimpose diagram C on diagram A. 


lll i 

1 ' 

•X’X*X’X*XvX-xvX": 

4 1 

X*X%*t*X<-X<*X'X*X*X 

5 1 


X*XvX*X*Iv.vXv.''.v. 

8 I 

1 i 

: 

1 , 1 

m 

111 

j 


141 

\15\ 

1 


18| 

191 



|22 

|B| 

1 

25 1 


i 

28 

129] 

; : 

1 1 

32 1 

Ilf] 

1 i 

: 

j 

36" 

37 1 

I 1 

j i 

[40 1 

• 

j: ; 

1 i 

42 1 

43) 

\ 

1 i 


|46| 

i47| 

1 1 

1 1 
1 1 

m 

51 

i 

i 

1 i 


54 1 

55 1 

1 

i 571 

1 1 
1 1 

1 

'' 

i&»»aoooowoooo6» 

1601 

1 i 

|6l| 

1 i 

1 I 

rill 


DIAGRAM B 
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. : 


i /T -1 

■ i 

1 

i i 

j : 

■ 

: j: 

: : i ^ ?■! 

i n ^ ^ 

155! 

i54i 


151 1 

50 II 1 

147 : 

|46| 


1431 

42 1 1 

1 40 1 


127 

|36i 


1 33 1 

32 1 ! 


|29l 

|28| 


1 1 25 y 

!| 23 I 

|22| 

1 1 

1 19! 

|l8l 

1 15 f 

i A :! 

1 14 1 


1 11 1 

lioll 

1 8 1 


5 I 

4 1 


1 1 1 


|[^[2 JQ |[61 1 66J[^6 JIJ 7 l[^ 

||171| 20lf21 |[ 43 II 42 11 24 
|^[ l i Ir H Il 37 l llil l 3 0 I I 3 1 |f33 
|[^2ir34l[35ir^r28 || 38 | | 39 | f25 

Ir4i' 


DIAGRAM C 


FINAL ENTRY 

DIAGRAM C 
SUPERIMPOSED ON 
DIAGRAM A 
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VI. An Algorithmic Perspective 

The arrow and bulging-out rules may be implemented to generate odd magic 
squares by an intermediate level computer scientist. Likewise the cross-diagonal 
elimination approach may be implemented. The focus would be on reducing the 
complexity and optimizing the execution in terms of time and memory. Pseudo code 
for the arrow rule may resemble the following: 

Start at Cell (1, («/2)+l), loop 1 to n 

Position loop value diagonally if possible using current row minus 1 and 
current column plus 1 as coordinates 

If not possible then arrow is considered blocked, then reposition to fall to 
theright using last row and current column plus 1 if the result falls within the 
square with coordinates (n,n) else position the successive integer in the cell 
directly below using row plus 1 and current column as coordinates 

VII. Open Questions 

We pose the question: Suppose n is an even number but it leaves a remainder of 2 
upon dividing by 4. Is it possible to have n X n magic square in this case? To a 
number theorist «=2(mod 4). 

VIII. Other Magic Figures 
Look at the magic triangle. 

1 



2 


10 


8 


3 
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Look at the magic column. 

1 4 5 

6 2 7 

8 9 3 

1 i i 


same sum 

Look at the magic star. 



IX. What kinds of questions do the magic triangle, column and star pose? Are some 
of these questions passed on to us through generations through folklore? 
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Abstract 

In this paper, we introduce and study essentially a class of multidimensional modified 
fractional calculus operators involving generalized Riemaim Zeta function and a general class 
of polynomials in the kernel. These operators are considered in the space of functions Mt (R "+ ). 
Some mapping properties and fractional differential formulas are obtained . Also images of 
some elementary and special functions are established. 

(Keywords : Fractional calculus operators/Mellin transform/Riemann Zeta function/General 
class of polynomials and multivariable H- function.) 

Introduction 

The H-fiinction of r-complex variables z\, defined by Srivastava and 

Panda' by means of a multiple Mellin Bames type contour integral is represented here 
in the following contracted form^ 


H{2u 


0, n : mi, «i; ;mr,nr 

,Zr]=H 

p, q : pu q\, ,Pr.qr 
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1 

(2710))' 


J 







( 1 ) 


where o = V- 1 

For the convergence, existence conditions and other details of the above 
multivariable H-fiinction, we refer to the book by Srivastava et. af'. 

Srivastava^ introduced and studied a general class of ijolynomials which is 
defined by 


k=0 ^5 


fMK 


A 


N=0,\,. 


( 2 ) 


The coefficients Aux W ^ > 0) being arbitrary constants, real or complex and M is an 
arbitrary positive integer. 

This general class of polynomials (2) unifies and extends a number of classical 
orthogonal polynomials such as Jacobi polynomials, Hermite polynomials, Laguerre 
polynomials, Gagenbauer polynomials, Bessel polynomials and several other classes 
of generalized hypergeometric polynomials. 

Throughout this paper we use some notations. As usual R and C represent the 
fields of real and complex numbers respectively. IC denotes the set of «-tuple real 
numbers, i?” non-negative real numbers and C" complex numbers. For brevity, we 

write for the product x^and x'’ for xf x^ with x =(xi, ,x„), X, = 

(Xi,... .,X„) e C” and p eC. Further we write i X] for the sum Xi + + X„. By (p+ we 

mean the positive part of a function (p defined by 


[0, if (p(x)<0. 


( 3 ) 


Goyal and Laddha'* introduced an extension to the generalized Riemann Zeta 
function which is represented here in the following slightly modified form. 


+ ^ 1, I z 1 < 1, Re(a) > 0 (4) 
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obviously when p. — 1 in (4) , it reduces to the following generalized Zeta function^. 

00 

(D(z,5,o)=X(«+«)‘^^", |z|< 1, Re(o) > 0 (5) 

M=0 


Modified fractional integrals : 

The multidimensional modified fractional integral operators and 

defined as follows 


r--n 


1 


= -Fjrn) &:zzs:'« 

Oa,(zi, s, c, g) 5^ [z {min (xj/U, ,x„/t„) - 1 dt (6) 

and fl.) = x,./,.)]; 


Oo,(z 2 , s, c, g) [z {1-max (xi/U, /!„)}“] /(O dt (7) 

for Re(p) > 0, where 

zi = z {min (xi/ii, ,x„/t„) - 1 \zi\< 1 (8) 

and Z 2 = z {1-max (xj/ti, ,x„lt„)}‘^, \zz\ < 1 (9) 

Special cases of the fractional integral operator (6) and (7) : 


The fractional integral operators (6) and (7) have a large number of special cases 
due to the presence of the general class of polynomials and generalized Riemann Zeta 
function Ob(z, s, c, g ) in the kernels of the integrals. We mention below a few of them 
for the sake of illustration. 

(i) If we set s = 0 and g = 0 in (6) and (7) then Ob(zi, s, c, g) reduces to unity 
and we get fractional integral operators recently studied by Goyal and 
Salim®. 
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(ii) If we set ^ = 0, g = 0, M = \, N = 0 and Ao,o = 1 in (6) and (7), then both 
general class of polynomials and Riemann Zeta function reduce to unity and 
we get the modified fractional integrals introduced by Brychkov et. af and 
studied by Tuan and Saigo* and Raina®. 

(iii) By expressing the general class of polynomials and Riemann Zeta function 
involved in (6) by their series form (2) and (4) and changing order of 
summation and integration, we get 


/+;« Zj t I • 

^=0 *=o k\g\ 


( 10 ) 


Tiji^cg+ka + l) 

r(p + l) * 


(iv) if we set « = 1 in (6) [ or in (10) ], we get 


(- H (g +s)'’z 

f=o i=0 ]c\g\ 


i+g 


( 11 ) 


r(|X -h eg + fa? H- l) j^+cg+ia^-n-cg-ka 


where is the well known Riemann-Liouville integral operator. 

(v) on setting o = 1 5 = 1, M = 1 and in (6), 

Riemann Zeta function and general class of polynomials reduce to hypergeometric 
function^ and Laguerre polynmial'“ respectively and we get 


I +;n 


a. 





ON MULTIDIMENSIONAL MODIFIED FRACTIONAL CALCULUS OPERATORS 23 1 


2 F 1 (1, c; 1+c; zi) Zl"' (z{min(xi//,, f(t) dt (12) 

fori?e(^i) > 0, |z,| < 1 

(vi) on setting o = \, s = -m {m = 1,2,3, ) in (6), generalized Riemann 

Zeta function is expressed in terms of Bernoulli’s polynomial^ and we get 




5” 


r()i + l)&, « 


J [min {xxitx, ,xjt„)-\l 




[z{min (xxItx, ,x„,/f„)}‘’]/(0 dt (13) 

for |log zx\ < 271 and B„+g+x(c) is Bernoulli’s polynomial. 

(vii) by dividing i?” for a fixed x e R" into « sub domains with zero- measure 
intersection 


R: = U {t e F; I (Xfc/fO ^ (V0)> 2=1 = j^k} 

k^l 


the multidimensional fractional integral operator can be expressed as a fimte 

sum of single integrals 




r(ii+i)i^xdx, 





( 15 ) 


Similarly by dividing i?” into sub domains with zero- measure intersection. 
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6 Rl I {Xklh)> ixjitj), j = 1, ,n ■, ji^k} (16) 

k=\ 


we obtain 


— 1 "5 

r{n+i)wac. 


X. ,^.<)]dl 

I 

(17 

Similarly, several (known or new) results may be obtained from (7). 

Modified fractional integrals of special functions : 

(i) Putting y(x) = x'^for e c" in (10) and making use of a known result ^ we get 

v,NM-.d,g d ^ y [^A^] (- (fj,)g (c + g) z'‘^^ 

h h k\g\ 


r{}J. + cg + ka + l)r(« -\i-cg-ka + 1£/|) 
r(p+i)r(«+|(^|) 


(18) 


provided Re{dj) > - I (J =1, ,«) and n + Re( |<il ) > Re( \i + eg + ka ) {k = 

(ii) Sety(x) = (ex" + bf H [ (ex" + 6) , , (ex" + b) in 

(10). We replace multivariable i^-function by its MelliurBames contour 
integral with the help of (1), collect powers of (qx^ + b) and apply the 
binomial expansion 

X! 




(19) 
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where {ax'"lb)<\ ; prg {cc^lb) < n 

Now, collecting the powers of:r and applying- a known result® we get 


I +;w 











g,6»=0 iir«0 A:!g!8!r(p + l) 


ff 


0,n' + 2:m,,n,; 

p + 2,<J + 2 . Py,(J\, — iPr^^r 


yPrxV-rb°r 


,(!,), (l -n + p-f eg + to- 6lv|-|4|;u 

(5 + X; a, , , a, ^ (2 - « - 5|vj; w, u. 







(20) 


(iii) 


Trm\n' 

Set / a:) = x“* 
known 




in (10), for h € R+ , where //(x) is well 


Fox's H-fiinction''. For details of this function one can refer to Srivastava et af. 

To evaluate the modified fractional integral of this function, we use the 
multidimensional Mellin inversion formula’: 

/W=(i/2«r (2i) 

for the Mellin transform 


Here, and in what follows, we use the notation 
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|•(T)+J00 j*t]+ic0 rX„+Z00 

i(x)-iQ0 * Jxfj-ico 


We shall require the following known results due to Raina^ : 

J x-/(min{x.'' , } )^& = (23) 



Re(sj) < 0 (j = 1, ,«) and 

|^x''7( max {x,\ a!" ))^ = — (24) 

•^1 


Re{sj) > 0 (/■ = 1, ,«), where /(/) denotes the one-dimensional Mellin transform of 

f[x). Now, in view of (21) and by application of (23) we can express 








°r{s)- 


-x’-'ds 


(25) 


where 


nrk+s^lMnrli-o.-^d#)} 

F(s)=-£i tL— 

nr{i-4,-Bjj|/A))nr{o,+4(H/4 

+I >n'+l 

Now operating on both sides and using a known result* we obtain 

Y:r-{A:-F;;[min(x‘, x;)]} 

_ V* (~ + g) '2 **T(p + eg + fe? + l) 

,=o 4=0 ^!g!r(p+i) 
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(- 1 )"'' ''’j:'" T{-\i-cg-ka- j^j) 

h{ 2 m)" (,^i. r(- s ) 



(26) 


Now , interpreting the R.H.S. of (26) as the iZ-fiinction by using (25), we get 


(xf , x„")]} 


« fAT/yl (~ + g)~' + eg + + 1) 

= l£ k\g\r{vi+\) 


x'^H 


jz+i.^r'+i 


lniin(x*, 





(27) 


provided that T;< 0 (/ =l,....,n), -lx| > Re(p + cg + to )> 0 (r- 0,1 ,[N/M \ ); and 


(i) A = I A, - +Y.Bj~ TBj > 0, or 

7=1 /=«'+! 7=1 7 = m'+l 


^ P ^ 


(ii) A=o, Re[i;»,-£i,l-o>'-9')/2+Mi2;^,-s;s, 

V;=i ;=i ; VJ=' >=i J 


> 1 ; 


Re{a^ < 1 - (Ixj A// 2 ) (/• = 1, ,n’), Re{bj) > -\x\BJh (j=\, ,m' \ 

Note that on setting s = 0 and g = 0 in equation (18) and (27), we get the result 
established by Goyal and Salim^. If we put s = 0, g = 0, ^ = 0, 7V^=0, M = 1 and 
A 0.0 = 1 in equations (18) and (27), we get the results established by Tuan and Saigo* 

and by Raina® respectively. Moreover , if.4j = 1 (/ = 1, ,p' ), Bj = 1 (/' = 1, ,q') 

and /2 = 1, then '(27) yields the corresponding formula for Meijer’s G-function 
obtained by Tuan and Saigo® 
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Modified fractional operators in space M, R " : 

Following*, let M^r" denote the space of functions /which are defined on R" , 
where t = (ti,...., x„) e R" . It is proved that / 6 M (/?"], if and only if, / can be 
represented as the inverse Mellin transform, 

(t)+J 0 O 

f(x) = (l/2ni)" jf'{s)x'‘‘ds (28) 

of a function /(s) infinitely differentiable and with compact suffort on ((x)-/oo, (t)+/oo). 

Theorem : Let /?e(p) > 0 ; X; + Re(dj) < I (j= I, n); Re(\i + eg + ka) + Re{ )+ |t| 

<n(k = 0, 1, , [N/M\ ){ord e c" and x € 7?", then the operator is a 

homeomorphism of the space onto itself. 

Moreover , it can be written in the form 


^ Igl (- iVL A, H (c + g)'‘ + eg + to + 1) 

‘=0 ^:!g-!r4i + l) 


1 ''* 1 " r(- p - eg - to - |rfj - 1 j| + n) 

( 2 ®) (t)-m r(-jc7| — |s| + Alj 


F’{s)x-’ds 


(29) 


Proof : Making use of (18), (28) and changing the order of summation and 
integration, we easily get (29). The interchange of order of integration is possible 
since/ ( 5 ) has a compact support . The function 


g. [^] ( Axj, (p)g (c + g )~' + eg + to + 1 ) r(n - p - eg - to - 1^1 - 15|) 

^=0 i=0 ^!^r!r(p+i)r(n-|£/|-|5|) 



ON MULTEDIMENSIONAL MODEFffiD FRACTIONAL CALCULUS OPERATORS 237 


has compact support and is infinitely differentiable on ((t) -i <x>, (t) + i oo) if so does 
/ (5). Hence belongs to Mx \ The continuity of the mapping 

/ /in Mx ) is obvious. 

Note that on setting g = 0 and j = 0 in (29), we get the results established by 
Goyal and Salim^ 
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Abstract 

A Bianchi Type K string dust cosmological model in General Relativity, is obtained. To get a 

determinate solution, we suppose s = X and assume a transformation di = odt where a is 

function of t and a is metric potential, e the rest energy density and X the string tension density. 

The physical and geometrical aspects of the model, are also discussed. 

(Keywords : cosmological model/relativity^ianchi type K /string) 

Introduction 

The choice for the study of Bianchi type IX cosmological model creates more 
interest because familiar solutions like Robertson-Walker universe, the de-Sitter 
universe, the Taub - NUT solutions etc. are of Bianchi type IX space-time. Bianchi 
type IX cosmological models include closed FRW models. 'These models allow not 
only expansion but also rotation and shear and in general are anisotropic. After the big 
bang, cosmic strings arise during the phase transition as the temperature goes down 
below some critical temperature’'^ 

'The general relativistic formalism of cosmic strings is due to Letelier'*'^. Stachel® 
has investigated massless strings. Baneqee et aV have investigated some 
cosmological solutions of massive strings for Bianchi type I space-time with and 
without magnetic field. Chakraborty and Nandy* have studied Letelier strings model 
for Bianchi type II, VIII and IX space-time. Chakraborty® has also investigated a class 
of cosmological solution of massive strings in Bianchi type IX space-time using a 
supplementary conditions a = a b" between metric potentials where a and b are 
functions of time-alone and a is a constant. Bali and Dave’® have investigated some 
special strings solutions for Bianchi type IX space-time using the condition e= ^ and 
a =e“' where a is metric potential. 

In this paper, we have investigated a general string dust solution for Bianchi 
type IX space-time using the transformation dx = a dt and the condition e = X, where 
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a is metric potential and function of /-alone, s is the rest energy density and X the 
string tension density. The physical and geometrical aspects of the model, are also 
discussed. 

We consider the Bianchi type IX metric in the form 

ds'^ = -dt^ + a^{t) dx^ + b'^it) dy^ + {h^ sin^ y + o^cos^y) dz^ 


-la^ cos ydxdz (1) 

The energy momentum tensor for string dust is taken as 
T/ = s v,.v'' - X x,.x-' 


with 


v'v. = -x'x, = -l (3) 

and 

v‘ Xi = 0 (4) 

where s = Sp + X, is the rest energy density for a cloud of strings with particles attached 
to them, 8p, the density of particle, X, the cloud strings tension density, v" the flow 
velocity vector and x" is the direction of strings. 

The Einstein field equation 

R!-\Rgi=-V (5) 


[using the unit in which 


StiG 


= 1] 


for the metric (1) leads to 
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^ai) P 1 

4b^ 


J b^ 1 3a^ 

b^ b'^^ b'^ 4 b* 


X 


a b 

ab 

1 o" 


+ — 

-j 

a b 

ab 

4 b* 


( 6 ) 

( 7 ) 

( 8 ) 


Eqn. (6) to (8) are three eqn. in four unknowns a, b, e and X. To get a determinate 
solution, we assume 

8 = X (9) 

Subtracting (7) from (6) and using the condition given by (9), leads to 


ab~ b'^ 4 b* 


( 10 ) 


Again subtracting (10) from (8), we have 


2- + ^ = 0 (11) 
b a 

where a dot (.) denotes di ff erentiation with respect to t. 

To obtain exact solutions of (10) and (1 1), we make use of the scale transformation 

dz = a dt (12) 

Using the transformation given by (12), eqn. (11) and (12) reduce to 



( 13 ) 
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222 :, + 4 = 0 

ab b a a 


where a dash denotes differentiation with respects to x. 
Eqn. (13) leads to 


dx 




W 


which on integration leads to 


h' = 



where a is the constant of integration. 
Eqn. (16) further leads to 


2b db 



f— T 


= dx 


using 




2-/a 


= 3 ' 


Eqn. (17) leads to 


23c/5 

2^/a3 



(14) 


(15) 


(16) 


(17) 


(18) 


( 19 ) 


which on integration leads to 
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( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 
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Using the transformation 




dx 


(29) 


in eqn. (28), we get 
d^A 


d^^ 

which leads to 


T =-A 


(30) 


sin"’ — = '¥ -2 tan"’ m 


where ^ , »j are arbitrary constants and 
'R = 2 tan’ (2aT + p) 

is obtained from (29). 

Now eqn. (30) leads to 

%o}lmx^ + 4aT(^+2p^ffl - £m^) + (jP - /w)(l+ P/w) 21 


(31) 


(32) 


A^a^- = 


(l + ) \4a\^ + 4apT + P^ + 1; 


(33) 


where i ,m axe arbitrary constants. 
Hence the metric (1) reduces to the form 

ds^ =- 


(l + m 

(4a^T^ + 4a px + p^ + 1 

1 

ia^imx^ +4ai:(j 

?+2p fm-.^»i’’) + (ip-/M)(l+pm) 2£ 


dx^ 


ia'^imx'^ +4ax (i 

[+2^£m-£m^] 

1 + (jP - /Mj 

I(l + p7«) 2i 

1 

l+ni^ 

) (4aV+4apT+p^+lj 

1 
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+ 


+ 


(at^ + Pt +y) dy^ + [(ax^ +Pt + y) sin^ y 

ia}imx^ +^ax{l + 2p im - i^nP' ) + (p - /«} (l + ^m)2i. 2 1 
(l+/w^)(4a^x^ +4apT+p^ +l) ^ 


dz^ 


2(8a^^/«T^4-4aT(^ + 2p^w-^m^) + (jp-ffl) (l + Pw) li 

(1 + (4a^T^ + 4apT + P^ +lj ^ ^ 


(34) 


Some Physical and Geometrical Features 

The scalar of expansion (9), the shear scalar (cs\ the spatial volume V and the 
density (s) are given by 


(8a^^m^+4aT(^ + 2p^/»-^/«^) + (P-/«)(1 +Pot)2^ 
( 1 + W ) (4a^T^ + 4apT + p^ +l) 


X 


2ax +P 

4a^mx+a 

_i < 

[l+2^m-m^] 


(ax^ + Px + y] 

4a^mx^ +2ax 


+ (P- 

w)(l + pm) 


4a^T + 2aP 
4 a^T^ + 4 apx + p^ +1 


(35) 


a 


2 


3 \b a j 


2 

3 


(8a^^m^+4aT(g4;2p^/K-^/«^) + (jp-/»)(l + pffl)2-g 
{l + /w )(4aV + 4apT + p^ +l) 


lax +P 


2(ax^ + Px+y) 


4 a^wT^ 


4 a^mt+a(l+ 2 pOT-ff?^) 

+2ax(l + 2p/«-m^)+(jP-»j) (l+p/n) 
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4ah + 2aP 
+4apT + +1 

V=b^a 


(36) 


= (at^+Px + y) 


(8a^^mT^+4at(^ + 2p^ffl-^/n^) + (p-w) (l + pw) 2^1^ 


(l + ff3^)(4aV + 4apT + P^+l) 


(37) 


dh 

ah 




lima. +/M^ + 1 
(l+»j^)(ax^ + Px +y) 


(38) 


Discussion 

If there are no particles attached to the string i.e. particle density ep = 0. Then a 
situation arises where energy density (e) is equal to the string tension density (A.) 
(from 6p = € - X,). This is the case of geometric string considered by Stachel®. This 
condition is physically relevant in cosmic explanation. 

The reality condition e > 0 is satisfied when 2 1 ma + + 1 > 0. There is a big 

bang in the model at x = 0. The scalar of expansion (0) is monotonically decreasing 
for X > 0 and 0 -» 0 when x -> oo. The model (34) represents shearing and non- 
rotating universe. Shear (o) is non-zero for 0 < x < oo. Since lim — ^tO. Hence the 
model (34) does not approach isotropy for large value of x. 

TT L lima+m^+l 

However, when x 0 then s = — r rx 

(l+TK^jy 

The spatial volume (V) is finite at x = 0 and it becomes infinite when x-> oo. e->0 
when x-» oo. Thus the model is essentially empty xmiverse when x-> oo. 
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Abstract 

The authors have considered a pair of dual series equations involving a pair of polynomials 
suggested by jacobi polynomials which are biorthogonal over the interval (-1, 1) and have 
obtained the unknown coefficients involved in the series equations. Their limiting case, leading 
to Konhauser biorthogonal polynomials have also been considered. 

(Keywords : Dual series equations/biorthogonal polynomials) 

Introduction 

In potential theory, many mixed boundary value problems reduce to the problem 
of solving of some dual series equations involving classical orthogonal polynomials as 
its kernels. 

In the same effort Srivastava and Panda* have considered the problem of solving, 
the following dual series equations involving Jacobi polynomials of different orders. 
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and 


I A 


(5 + « + .^ + l) 


(a + « + £ + 1) 


P 


fl-— 1 
0 ) 
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Where, c>0, f is an arbitrary non-negative integer, fif) and g(x) are prescribed 
functions. 
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and 


/,={x : 0 <x <>'}, /2 = {x : < X < c}, (3) 

and in general, 

Min {a, p, y, 5, X, [i, £ , a} >-l (4) 

For this purpose Srivastava and Panda' used the technique, of Noble^ with some 
modifications. The dual series equations (1) and (2) happened to be generalizations of 
those considered earlier by Lowndes^, Noble^ and Dwivedi & Trivedi"*. 

Recently, Madhekar and Thakre* considered the biorthogonal system of 
polynomials {/^^“’^^(x; A:)} suggested by jacobi polynomials over the 

interval (-1, 1) with weight fimctionH'(x) = (l-x)“. (1+ x)^ 

We know that 


(x ; 1) = (x;l) = Pi“-»(x) (5) 

Thus, in view of above, in the present paper the authors have considered the 
problem of determining the sequence {A„} satisfying the following dual series 
equations : 
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1 ;k 
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= g(x),Vx6l, 


( 7 ) 


Where c > 0, £ is an arbitrary non-negative integer, / (x) and g (x) are prescribed 
functions, and in general, 
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Min {a, p, y, 5, X, ja, p, a} >-l 

The method employed to solve (6) and (7) shall be analogous to that of used by 
Srivastava and Panda' with necessary modifications. 

In view of (5), the dual series equations (6) and (7) are generalizations of (1) and 
(2), considered by Srivastava and Panda', and hence are generalizations of those 
considered earlier by Lowndes^, Dwivedi and Trivedi"*. 

It is interesting to note that : 



( 8 ) 


(9) 


( 10 ) 


( 11 ) 


( 12 ) 
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In view of (8) and (9), by suitably appealing to the principle of confluence, our 
results in this paper would readily yield the corresponding results for dual equations 
involving series of Konhauser polynomials. 

Preliminaries 

We will use the following results which may be proved easily by use of 
differential recurrence relations : 

(i) Fractional integrals : 
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(ii) Derivative formulae : 
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|(P+/M+w+l) r 



1 ,k ,m>0 

c ; 


( 16 ) 


(iii) The biorthogonality property in slight modified form : 

f 1 --T 1 - 1 - — 

ovc; 

(1 + a + Icn) (1 + p + «) 

= —= p= 5^,a=P>-l (17) 

«! (l + a + p + «)(l + a + p + « + ^) 

Where 5^1 is liie Kronecker delta. The above results can be proved easily by use of 
above results and following the method of Madhekar and Thakre^. 

Corollary of (17) : 



0<>/^c 
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Where y is arbitrary. 

proof : To prove (18), for m 9^: n we have, from (17) 

|%“j^l - - y ;)S: = 0 (19) 

Similarly 

C 7 rY 7 9r ^ f Jv \ 

f x“ 1-- l-—;k 1 ;k dx=0 (20) 

^0 { c) \ c J K ^ J 

For 0 <y < c we add (19) and (20) and write it in following form : 




( 22 ) 

y \ cj 

So that (21) can be written as : 

Clim, n; y) + n (n, m; y) = 0 m^n 
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or, 


Q.{m, n, y) = - Q. {n, m, y) = Q m^n (23) 

From (22), we notice that Q (m, n;0)= Q (m, n;c) = 0 

For m^n and (23) holds for every arbitrary y, so we conclude that £2 {m , «; y) = 0, for 
ally e [0, c]. 

This proves (18). 

In view of (5) it is further to note that for ^1. (18) reduces to the modified 
orthogonality relation for jacobi polynomials, given by Srivastava and Panda’, (18) 
shall be most useful property in our investigations ; 

Multiplying-Factor Technique 

I 

We multiply (6) by x(<|)-x)'"'^^"’ where p is an arbitrary constant and w is a 
suitable non-negative integer, and integrate both sides with respect to x over (0, <{)), 
then by use of integral formula (13), we obtain : 

\y +A:(« + -g)H-l]|[l-Hx + ^(« + ^)]| + ^ 

"-0 " [p + k(n + ^) + 1] l+a + m + p + k(n + tf] 


. j = J - xr^^f{x)clx 

Now Difif. (24) both sides m times w.r.t. (j), we get by use of (15) : 
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Where 0 <({)<>', a >-l and»2+;?>0 

Now when we multiply (7) by 1 - - and differentiate the resulting equation / 

I cj 

times with respect to x, then using formula (16), we obtain ; 
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Where / is a non-negative integer. 

Now multiplying (26) by (x-(|))*^''' and integrating both sides with respect to x 
over the interval ((j), c) and using the formula (14) we obtain : 
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Wherey < (j) < c, p - / > -1 and /z + / > 0. 

Now if /) and /j are so choosen that the coefficients of and 

become identical and these polynomials are of the same order, then by applying (17) 
and (19), unknown coefficients {A„} can be determined. Below we discuss all such 
situations of interest : 
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The special Case p = a and a = p 

In the special case p = a and or = p, we get, p = y - a and h = 8 - \i, the 
equations (25) and (27) reduce to the form : 
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provided that 

a+p=y+5=X+p 

Now appealing to (17), (18) and (30), we get our result contained in : 

Theorem : For c> 0, let the sequence {A„} be defined by the dual series equations 
(6) and (7). 

Then for integers i, m, 8>0 
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(X-Y + /)V S + 




(31) 


Where for convenience, 


F(,j)) = d; I j V (<}) - 


and 


<p 



g{x)dx 


(32) 


(33) 


Provided that (30) holds and 

a + p + l>-l, y + » 2 >a >-l, 5 > |r - / > 1 (34) 

This theorem reduces to that of Srivastava and Panda* for A: = 1 and hence also 
generalises to that corresponding results of Dwivedi and Trivedi'*. 
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Abstract 

Results of an experimental study, earned out on radiations from neuronic shaped planar 
microstrip patch and array antennas are presented in graphical from. It is concluded from this 
study that by applying scale-modelling principle on such radiators, radiation properties of 
human brain neurons can be estimated. 


(Keywords : microstrip antenna/neurons/radiations from brain/ scale-modelling principle) 


Introduction 

The radiation properties of microstrip antennas have been extensively studied in 
recent years due to their practical advantages over conventional antennas including 
their lightweight, ability to conform to the host object and relatively easy and 
inexpensive fabrication compared to other antennas^"^. Earlier work related with 
radiation properties of microstrip antennas was limited to microstrip patches of regular 
shapes particularly rectangular, circular, triangular or elliptical geometries. In recent 
years, interest of scientific community has started shifting towards the analysis of 



260 


DEEPAK BHATNAGAR et. al 


planar arbitrary shaped microstrip antennas^^ by applying different mathematical 
approaches. 

Human body in general contains more than 10^^ neurons and most of them are 
concentrated in human brain and eyes^’^. The sizes of these neurons vary from 0.02 
micrometer to 100 micrometer or more. The shape, size, dimensions and properties of 
neurons differ from one portion of the brain to other. Most of the neurons are made of 
conducting materials. Due to continuous variation in the concentration of sodium, 
potassium, calcium and chlorine ions in brain, positive and negative voltages of 
intensity varying from -70 milli to approximately +60 milli volts are constantly 
generated^. It is worth mentioning here that in recent times, these voltages have been 
amplified thousands of time and are used in a device known as "MIND SWITCH" 
which is proving useful for handicapped and severely disabled persons. With the help 
of this device, these people can now operate electronic and electrical devices by using 
only brain signals without moving any part of his/her body^^'^^ When these 
alternating positive and negative voltages are applied to conducting neurons, they 
generate electromagnetic waves whose frequency depends upon the shape and other 
parameters of the neurons. In the same way, these neurons receive electromagnetic 
waves of appropriate frequency if available in the surrounding. Some AC signals of 
brain frequency spectrum have been measured by available instruments and found 
them in the range of few hertz. The voltages produced by these signals have been used 
to obtain Electroencephalography (EFG) signals from the brain^^. 

The concept of designing of neuronic patch microstrip antennas (NPMA) to 
analyse the radiation and reception behaviour of human brain neurons was developed 
from two dimensional photographs of human brain neurons. The shape of a neuron of 
human hippocampus as seen by Barinaga*^, is similar to an irregular patch antenna. 
Human neurons have three-dimensional geometry but at this stage of work , to 
simplify our problem we have considered their cross-section only and developed them 
as microstrip patch antennas having small thickness. By applying scale modelling 
technique^"^'^^, dimensions of the cross section of a human brain neuron are converted 
into an antenna structure capable to operate in C band of electromagnetic frequency 
spectrum. However the thickness of neuron is not scale modeled due to non- 
availability of dielectric substrates of different thickness in local market, which are 
required for fabrication of microstrip antennas. This is perhaps the first reported work 
on NPMA structure hence scope for incorporation of third dimension Le. thickness of 
neurons exist for future work to get more realistic picture regarding radiation 
performance of human neurons. This study on the radiation properties of neuronic 
antennas can be proved helpful in improving the performance of devices like "Mind 
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switch" etc. by indicating the direction of maximum intensity and nulls of the 
radiation patterns of brain radiators. 

Design Details and Results 

There are several shapes of neurons found in human brain'®. The dimensions of 
human brain neurons are so small that an antenna structure with exact dimensions of 
an actual neuron cannot be fabricated for experimentation. Considering this fact, scale 
modelling principle'"'''® is applied and dimensions of cross section of a human neuron 
are increased so that neurons may be visualised as neuronic microstrip antennas with 
resonant frequency lying in C- band (4 to 8 GHz). The geometry with coordinate 
system of one such neuronic microstrip antenna developed in our laboratory is shown 
in Fig. 1(a). The patch is located on one side of dielectric substrate (glass epoxy) while 
a very thin copper layer lies on other side of the substrate as ground plane. The 
substrate thickness is considered to be h and relative dielectric constant s^. The 
thickness of neuron (along z axis) is not scale modelled in the present paper. In fact, 
human body has hundreds or thousands of neurons associated in complicated three- 
dimensional networks. Instead of preparing a network at this stage of work, we 
initially prepared and observed tire behaviour of simple two-dimensional linear and 
planar array arrangements of neuronic patch antennas. A. two elements array of 
neuronic patch microstrip anteimas as shown in Fig. 1(b) and a 2x2 elements planar 
array of similar antennas are also prepared for the present study. Both patch as w’ell as 
array antennas are designed on a commercially available glass epoxy substrate (Sr = 
3.55) having 4' x 4' ground plane. 

Several approaches are available to analyse arbitrary shape structure In the 
absence of regular boundaries, most of the techniques in general, are either based on 
rigorous mathematical treatments or they require sophisticated simulation softw'are. 
Theoretical investigations of regular shape microstrip antennas require equivalent 
outward extension to incorporate the effect of fringing fields at the peripheries of the 
patch'*. The extension depends upon the planar dimensions of the patch, substrate 
thickness, substrate relative permittivity and field distribution at its peripheries. It is 
difficult to estimate the extension at the peripheries of an arbitrary shaped patch. 
Palamsamy and Garg'^ proposed that for a given irregular shape, extension parameter 
available for a shape, closest to that of given shape could be used. The extension 
parameter involved with a rectangular geometry is well known'*. Therefore a 
rectangular geometry OABC as shown in Fig. 1(a) just fitting outside the designed 
neuronic patch is considered to incorporate extension in different directions. 
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The length L of this rectangle OABC is 2.0 cm while width fF is 1.2 cm. The 
substrate thickness h and substrate relative permittivity Zr, applied for the present 
analysis are 0.16 cm and 3.55 respectively. The computed extension parameter AZ, is 
0.469cm. By applying a computer programme named 'DEEREK. CPP' prepared in 
Window-95 based programming language, analysis of this rectangular geometry is 
carried out by applying Model Expansion Technique'®. Any point on the periphery of 
rectangular patch can be used to feed it if a suitable (A, / 4) transformer is attached 
between that point and feed line. Therefore point 'P', shown in Fig. 1(a), lying on the 
periphery of rectangle OABC may be used for feeding this rectangular structure. The 
computed input impedance (ZJ of this rectangular antenna structure at point 'P' is 
(1 19.1 + j 2.87) ohm while its computed radiation efficiency is 79.1 1%. An argument 
similar to that of Palamsamy and Garg'^, regarding extension parameter of arbitrary 
shape structures, can be also be made here for searching the feed location on a 
neuronic patch antenna. The argument used in this paper is that for feeding a given 
irregular shape, feed location available for a shape, closest to that of given shape 
could be used. It can be seen from Fig. 1(a) that rectangle OABC just touches the 
neuronic patch at the point 'P' (1.10 cm, 0). Therefore point 'P' is used for feeding 
NPMA structure. For this purpose, computed input impedance of equivalent 
rectangular patch anterma (Zin) at point 'P' is equated equd to that of a NPMA 
structure. To match this input indepdance with that of feed line (50 ohm), a (7, / 4) 
transformer is designed and attached between point 'P' and the feed line. The 
designing of this transformer and computation of width of feed lines is carried out by 
applying another computer programme 'DEEOSA. CPP' prepared by using the 
relations of Seff and Zc discussed by Bahl and Bhartiya'*. Feed arrangements for two 
elements array structure and 2x2 elements planner array antenna are designed 
following the same procedure. For two elements array structure, feed arrangements is 
shown is Fig. 1(b). All the structures are fed through SMA coaxial connectors and 50- 
ohm coaxial line. The measmed E’-plane ((j) = 0°) radiation pattern of NPMA structure 
is compared with that of its equivalent rectangular microstrip anterma (RMA) 
structure Fig. 2. These two patterns are nearly alike in nature. Both structures have 
maximum radiation intensity in end-fire direction ((j) = O'*). The measured 3dB beam 
width of NPMA structure is 70° as shown in Fig. 2, which is significantly smaller than 
that of RMA structure (110°) computed under identical conditions. 

The radiation patterns, return loss, resonance frequency and input impedance of 
both neuronic patch as well as its arrays are measured at "Communication Systems 
Group, ISRO Satellite Centre, Bangalore" and "CARE IIT, New Delhi". The return 
loss, resonance frequency and input impedance measurements are carried out by using 
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-Measured NPMA Patterns 



Fig. 2-Comparison between measured £-plane patterns of NPMA structure and 
computed £-plane pattern of its equivalent RMA structure. 

a Vector Network Analyser and its associated computer programmes. During 
impedance measurements, due consideration was given to the accuracy enhancement 
techniques^® to correct effective directivity, effective source match and frequency 
tracking errors. The measured resonance frequency of patch antenna using network 
analyser is around 6.68 GHz, which is little higher than that of its equivalent 
rectangular microstrip antenna (5.89 GHz). The measured return loss for a neuronic 
patch antenna is (-12.48) dB at 6.68 GHz as shown in Fig. 3(a) while the computed 
return loss of its equivalent rectangular microstrip antenna is (-17.77) dB at 5.89 
GHz. The variation of return loss with frequency for two elements array and 2x2 
elements planar array are shown in Fig. 3(b) and 3(c) respectively. 
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Fig. 3(a)-Measured return loss of a neuronic patch antenna. 
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Fig. 3(b)-Measured return loss of two elements array of neuronic patch antennas. 
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Fig. 3(c)-Measured return loss of 2 x 2 planar array of neuronic patch antennas. 
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Fig. 5(b) 

Fig. 5(a)— Radiation patterns of neuronic patch antenna with different (ji planes. 

Fig. 5{b)-£-plane radiation patterns of a two-elements array of neuronic patch antennas. 
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The measured resonance frequency for a two elements array is around 6.725 GHz 
while for four elements 2x2 planar array it is around 6.935 GHz. The measured input 
impedance for a two elements array as a function of frequency is shown in Fig. 4, 
which indicates a significant mismatch between radiating elements and feed network. 
The mismatch between radiating elements and feed network as well as difference in 
measured and computed resonance frequency and return loss values are on the 
expected lines since the feed network for NPMA structures is designed by considering 
smooth boundaries of regular rectangular patch antenna but in reality, the considered 
geometry of NPMA radiator bears an irregular shape though it is close to a rectangular 
shape. 

The radiation patterns of NPMA structure and its two elements array are shown in 
Fig. 5(a) and 5(b) respectively. These patterns were measured by placing transmitting 
test antenna (NPMA) and fixed receiving horn antenna inside an anechoic chamber. 
The test antenna was mounted on an arrangement lying on a one-meter diameter 
circular wooden platform at ISRO Centre, Bangalore. The mount holding test antenna 
was capable in rotating in full range of colatitude angle (6) varying from 0° to 360° 
and in different azimuthal ((|)) planes. The separation between transmitting and 
receiving antennas was kept around two meters. With the help of a sweep generator, 
6.68 GHz frequency signal was applied and radiation patterns for patch antenna were 
measured by changing angle 0 from 0° to 360°. (|) = 0 and 90° pattern are sufficient to 
understand the radiation behaviour of NPMA structure. However radiation patterns 
plotted on polar graph are also reviewed in ({) = 45° plane as shown in Fig. 5(a) to get 
three-dimensional feeling. The distribution of radiation intensity in upper hemisphere 
is quite similar for different ^ planes. In (j) = 0 and 45° planes, maximum radiation 
intensity is in broadside direction while in (j) = 90°, maximum radiation intensity is in 
9 = 45° direction which are about 3dB more than that in 0 = 0° direction. The 
radiations in lower half are due to surface waves generated by antenna, which are 
undesired radiations and needs to be suppressed. The radiation intensity in lower half 
is much smaller than that radiated in upper half. For a two elements array antenna, E 
plane radiation pattern is measured experimentally in (j) = 0° plane at 6.725 GHz and is 
shown in Fig. 5(b). The main lobe in intensity distribution curve is again in broadside 
directions followed by side lobes. The side lobe level is about 15 dB down in 
comparison to the main lobe. Radiation intensity in lower half is again smaller than 
that radiated in upper half. 

Discussion and Conclusion 

Kirkup and co-workers" for the first time amplified the electromagnetic waves 
generated by human brain neurons in 1997 and found some interesting results in the 
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present communication, instead of observing radiations by actual human neurons, we 
have scale modelled the dimensions of a model of human brain neuron and converted 
it into a planar microstrip antenna operating in C band of frequency spectrum. Latter 
some planar arrays are designed and their radiation properties are measured 
experimentally. The expected radiation behaviour of human brain neurons is tried to 
explain with the help of measured data, of neuronic microstrip patch and array 
antennas.' 

If we decrease the size of the modelled neuronic antenna to that of actual human 
brain neurons (0.02 micron to 1 00 microns), its effective frequency comes in the range 
of lO''* Hz. This is the frequency band of electromagnetic waves radiated by the Sun 
and other heavenly bodies and received on the earth after filtering through the 
electromagnetic windows around earth's atmosphere. The neurons of eyes receive 
these electromagnetic waves and convert them into respective colours using relevant 
portions of brain. In the same way, neurons of different parts of the brain also receive 
and generate electromagnetic waves whose frequencies correspond to the shape, size, 
dimension and other properties of the neurons. 

The measured ^-plane radiation patterns of a single neuronic patch microstrip 
antenna are drawn in Fig. 5a. These patterns are similar to those of other regular 
shaped microstrip antennas. Different (j) planes are considered to visualize three- 
dimensional radiation patterns. It is observed that patterns are nearly uniform in all the 
directions. One can expect that in the same way human neurons might be radiating 
uniformly in all the directions. No back or side lobes are observed in measured 
radiation patterns of NPMA structure. Due to design considerations, NPMA structure 
can radiate only in upper hemisphere while an actual neuron will radiate in all the 
directions. The 3 dB beam width of neuronic patch antenna is little less than that of 
corresponding rectangular patch antenna operating under similar conditions. The 
measured E plane radiation pattern ((}» = 0°) of two elements neuronic patch array 
antenna excited in same phase is shown in Fig. 5b. The pattern has a well-defined 
main lobe followed by side lobes on both sides of main lobe. These side lobes are 12 
to 14 dB down to the main lobe. The present model is a simplified model still one can 
have a feeling that human neurons arranged in complex three dimensional network 
might be radiating/receiving electromagnetic signals in/from the surroundings in a 
similar way. 

The resonance frequency of a neuronic patch antenna is found 6.68GHz. It was 
expected that it would be significantly close to (5.89 GHz) i.e. of considered 
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rectangular patch aotenna. The expectation was again based on the fact that geometry 
of designed neuronic antenna structure was nearly fitting into the considered 
rectangular patch geometry. However a significant deviation of 12% between 
observed and computed resonance frequency is recorded. In the same way, input 
impedance and return loss results are showing a significant mismatch between 
radiating patches and feed line. The reason behind this difference is that input 
impedance of NPMA structure at point ’P' is equated equal to that of its equivalent 
RMA structure and feed lines are designed accordingly. This arrangement is not found 
fully suitable for feed point selection. For proper matching between antenna and feed 
line, actual input impedance at point 'F should be computed accurately by applying a 
simulation software and feed lines must be designed accordingly. Possibilities of error 
due to manufacturing limitations during designing of antenna structures cannot be 
ruled out. 

As mentioned earlier, the signals generated by human brain neurons can be 
amplified several thousand times by means of a device now known as Mind Switch^®. 
Research work in this direction is progressing^ The work will open possibilities of 
directing electromagnetic waves of particular phase and frequency towards a special 
part of brain and hence in modifying the behaviour of human being as brain neurons 
controls several characteristics of a person. 

This is the first attempt, where expected radiation performances of a human brain 
neuron are observed by modelling it as planar microstrip patch and array structures 
operating at microwave frequencies. By applying available simulation software to 
locate correct feed point on antenna geometry and by designing several such radiating 
structures in the form of three-dimensional networks, more realistic results may be 
obtained. 
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